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ABSTRACT 

A  transverse  resonance  mode  matching  technique  has  been  developed  to  analyse  passive 
microwave  and  millimetre-wave  waveguides  and  components.  This  technique  possesses 
superior  computational  efficiency  when  compared  to  more  general  approaches  such  as  the 
finite  element  method.  This  advantage  is  particularly  useful  for  analysing  broadband 
components  used  in  electronic  warfare  systems.  In  this  report,  the  theory  behind  this 
transverse  resonance  analysis  is  presented  in  detail.  Theoretical  results  for  several 
waveguiding  structures  are  presented  and  compared  with  experimental  or  published  data  to 
verify  the  analysis. 


RELEASE  LIMITATION 


Approved  for  Public  Release 


Department  of  Defence 


Aecsesioa  for 

KflS  QM&I 
DIIC 

Un  a.nn  o ;  3jt  c  e  d 
Just  if 1 cat  if 

□ 

□ 

n  . 

■ 

By - - - - 

Bistrlbctlc.n/ 

Avsilfcbllity  Cosoa 

Sist 

1'' 

Avail 

Spaa 

lal  1 

t 

DEFENCE  science  AND  TECHNOLOGY  ORGANISATION 


DSTO-RR-0027 


Published  by 

DSTO  Electronics  and  Surveillance  Research  Laboratory 
PO  Box  1500 

Salisbury  South  Australia  5108 

Telephone:  (08)  259  5181 
Fax:  (08)  259  5938 
©  Commonwealth  of  Australia  1995 
AR-009-211 
February  1995 


APPROVED  FOR  PUBLIC  RELEASE 


ii 


DSTO-RR-0027 


Transverse  Resonance  Analysis  Technique  for 
Microwave  and  Millimetre-Wave  Circuits 


EXECUTIVE  SUMMARY 

An  efficient  and  accurate  method  for  analysis  of  broad  bandwidth  microwave  and  millimetre-wave 
devices  has  been  developed.  This  report  describes  the  method  in  detail  and  the  results  obtained  by 
applying  the  method  to  typical  microwave  structures  are  compared  with  experimental  or  published 
data  to  verify  the  accuracy  of  the  method. 
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1  INTRODUCTION 


1.1  Purpose  of  this  study 

Computer-aided  design  software  for  microwave  circuits  is  now  readily  available,  at 
considerable  cost,  from  companies  such  as  Hewlett  Packard,  Compact  Software,  and  many  others. 
However,  these  software  packages  mostly  are  directed  towards  analysing  planar  integrated 
circuits  using  microstrip,  and  variants  such  as  coplanar  waveguide  but  do  not  adequately  analyse 
other  structures  in  use  especially  at  mm-wave  frequencies.  Structures  such  as  rectangular  or 
circular  waveguide,  finline,  dielectric  waveguides  are  only  analysed  in  an  approximate  manner 
(many  programs  are  dominant-mode  only)  if  at  all.  More  specifically,  some  of  the  limitations 
have  been  identified  to  include:  [1] 

1.  the  inability  to  predict  the  onset  of,  and  the  consequences  of,  higher  order  propagating 
modes. 

2.  the  lack  of  design  information  for  popularly  emerging  'new'  transmission  media,  e.g. 
finlines. 

3.  the  lack  of  design  information  concerning  active  devices  at  millimetre  wavelengths. 
(Low  frequency  design  models  cannot  be  successfully  extrapolated  because  of  the 
increasing  importance  of  parasitic  elements  and  the  distributed  nature  of  the  problem.) 

4.  the  lack  of  design  information  relating  to  discontinuities  in  transmission  media  at 
millimetre  wavelengths,  and 

5.  the  inability  to  predict  the  interaction  of  circuit  elements  which  are  in  close  proximity 
to  one  another. 

The  main  areas  of  weakness  of  existing  commercial  software  are  thus  in  the  areas  of  conventional 
waveguides,  finline,  dielectric  waveguides,  transitions  between  various  structures,  and  the 
capability  to  explore  new  waveguiding  structures.  These  problems  are  slowly  being  overcome 
through  the  use  of  field-theory  based  simulators  [2]  as  a  replacement  for,  or  in  conjunction  with, 
circuit-theory  based  simulators.  However,  even  these  simulators  are  mostly  oriented  toward 
planar  circuits  because  this  is  where  the  greatest  commercial  demand  lies.  The  exception  is  some 
software  (mostly  based  on  the  finite-element  method)  such  the  High-Frequency  Structure 
Simulator  from  Hewlett  Packard.  While  this  is  sufficiently  general  to  analyse  many  circuits,  it 
requires  a  powerful  computer  to  analyse  even  relatively  simple  circuits.  There  is  always  a  trade 
off  between  methods  that  are  relatively  simple  to  implement  but  require  a  lot  of  computation 
power  and  methods  that  are  highly  efficient  but  either  lack  generality  or  are  complex  to 
implement. 

This  report  describes  an  approach  which,  while  more  efficient  than  general  field-theory  based 
methods  such  as  the  finite  element  method,  still  enables  a  wide  variety  of  circuits  to  be 
analysed.  The  philosophy  has  been  to  pursue  the  middle  ground  -  that  is,  a  method  which  is 
general  enough  to  solve  a  variety  of  practical  problems  that  occur  particularly  in  millimetre- 
wave  circuit  design,  but  at  the  same  time  is  efficient  enough  to  be  used  for  routine  design. 

The  prime  motivation  for  this  work  has  been  that  many  electronic  warfare  applications  require 
extremely  broadband  components  (multi-octave)  and  the  constraints  in  the  design  of  these 
components  are  much  more  severe  than  for  narrow  band  components.  For  example,  in  the  design  of 
a  broadband  amplifier  or  mixer,  there  is  a  direct  trade  off  between  gain  (or  conversion  loss)  and 
bandwidth.  Moreover,  high  performance  for  e.g.  direction-finding  may  require  a  high  degree  of 
gain  flatness  and  matching  between  individual  components.  This  combination  of  broad 
bandwidths  and  high  precision  can  only  be  economically  achieved  with  accurate  computer 
analysis  techniques  which  are  specifically  designed  for  these  tasks. 
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1.2  Review 

A  wide  variety  of  numerical  techniques  for  the  analysis  of  microwave  and  millimetre-wave 
structures  have  been  developed  over  the  years  and,  aided  by  the  rapid  advances  in  computing 
technology,  many  are  beginning  to  now  reach  a  level  of  maturity.  However,  no  single  technique  is 
suited  to  all  problems  and  this  situation  is  unlikely  to  change.  At  present,  one  must  choose  a 
technique  (or  software  package)  that  is  appropriate  to  the  problem  being  solved,  just  as  a 
tradesperson  chooses  the  suitable  tool  based  on  training  and  experience.  What  will  happen  in 
the  future  is  that  design  software  will  have  a  range  of  techniques  built  into  the  one  package  and 
the  user  will  be  largely  unaware  of  which  technique  is  being  used  because  it  will  be  the  software 
rather  than  the  user  that  will  choose  the  most  appropriate  technique. 

1.2.1  Analysis  techniques 

Analysis  techniques  have  been  described  in  detail  in  recent  books  edited  by  Itoh  [3], 
Sorrentino  [4],  and  Yamashita  [5].  These  methods  are  summarised  in  Table  1.1  according  to 
the  chapter  headings  in  the  books  and  grouped  to  show  where  the  methods  are  common  to 
more  than  one  book.  There  is  also  a  certain  amoimt  of  overlap  between  many  techniques  e.g. 
the  transverse  resonance  technique  uses  mode-matching  which  in  turn  is  fundamentally  a 
moment  method.  Methods  other  than  transverse  resonance  and  mode  matching  will  not  be 
discussed  further  in  this  report. 


Table  1.1  Summary  of  analysis  techniques. 


Itoh  [3] 

Sorrentino  [4] 

Yamashita  [5] 

The  Finite  Element  Method 

Finite-Element  Method 
Finite-Difference  Method 

The  Finite-Element  Method 

Boundary-Element  Method 

The  Boundary-Element  Method 

The  Transmission  Line  Matrix 
(TLM)  Method 

Transmission-Line  Matrix 
Method 

Planar  Circuit  Analysis 

Planar-Circuit  Approach 

Spectral  Domain  Approach 

Spectral-Domain  Approach 

The  Spectral  Domain  Method 

The  Method  of  Lines 

Method  of  Lines 

Integral  Equation  Technique 

Method  of  Moments 

The  Mode-Matching  Method 

Mode-Matching  and  Field- 
Matching  Techniques 

The  Mode-Matching  Method 

The  Point-Matching  Method 

Transverse  Resonance  Technique 

Transverse  Resonance 

Techniques 

The  Waveguide  Model  for  the 
Analysis  of  Microstrip 
Discontinuities 

Generalised  Scattering  Matrix 
Technique 

The  Wiener-Hopf  and  Modified 
Residue  Calculus  Techniques 

The  Geometrical  Theory  of 
Diffraction 

The  Equivalent  Source  Method 

Asymptotic  Expansion  Methods 

The  Beam  Propagation  Method 

2 


DSTO-RR-0027 


1.2.2  Origins  and  background  of  current  method 

A  brief  summary  of  work  that  has  influenced  the  current  development  follows.  This  is  not 
an  exhaustive  review  but  summarises  key  developments.  The  beginnings  of  the  method  can 
be  found  in  the  early  work  on  microwave  network  theory  of  Montgomery,  Dicke  and 
Purcell  [6],  and  Marcuvitz  [7]  in  the  well  known  MIT  Radiation  Laboratory  Series  and  also 
a  little  later  in  Altschuler  and  Goldstone  [8]  who  developed  network  representations  of 
waveguide  obstacles.  Clarricoats  and  Oliner  [9]  derived  a  transverse-network 
representation  for  hybrid  modes  in  inhomogeneously  filled  waveguides  using  a  method 
employed  previously  in  the  representation  of  slotted  waveguides.  Felsen  and 
Marcuvitz  [10]  described  a  generalised  derivation  of  the  transverse  field  equations  and 
modal  representations  of  the  electromagnetic  field  in  transversely  inhomogeneous  regions. 
Kerns  [11]  discussed  concepts  and  conditions  underlying  the  establishment  and  use  of 
immittance-  and  scattering-matrix  descriptions  of  waveguide  n-ports.  Peng  and  Oliner  [12] 
discussed  a  class  of  open  dielectric  waveguides  which  is  of  direct  importance  to  the  areas 
of  integrated  optics  and  millimetre-wave  integrated  circuits  and  presented  a 
mathematical  formulation  based  on  a  rigorous  mode-matching  procedure.  Bornemann  and 
Arndt  [13]  used  transverse  resonance  to  calculate  the  characteristic  impedance  of  finlines 
with  up  to  three  slots  by  a  rigorous  hybrid-mode  analysis.  Masterman  and  Clarricoats  [14] 
described  a  computational  method  for  solving  a  wide  range  of  transverse  waveguide 
discontinuity  problems  and  showed  that  in  some  cases,  the  solution  is  found  to  be  sensitive 
to  the  way  in  which  infinite  series  of  field  functions  is  truncated.  They  further  showed 
how  the  optimum  form  of  truncation  can  be  determined  for  many  configurations  of  practical 
importance.  Ping  and  Jingfeng  [15]  investigated  the  filter  characteristic  of  NRD 
waveguide  by  combining  a  network  approach  with  mode  matching  theory  in  an 
application  of  the  method  of  Peng  and  Oliner.  Bates  and  Ko  [16]  used  multi-modal 
admittance  matrices  to  represent  the  coupled  radial  regions  in  a  waveguide  diode 
mounting  structure. 

1.3  Scope  of  this  study 

The  analysis  described  in  this  report  is  a  full  three-dimensional  mode  matching  analysis  using 
generalised  transverse  resonance  to  analyse  the  constituent  two-dimensional  cross-sections  in 
passive  waveguide  structures.  The  structure  is  assumed  to  be  lossless  and  consisting  of 
homogeneous  and  isotropic  cuboidal  sub-sections  which  may  be  dielectric  filled. 

1.4  Overview  of  contents 

Chapter  2  describes  the  technique  in  detail,  while  chapter  3  gives  results,  including  dispersion 
plots  for  a  range  of  two-dimensional  structures  and  S-parameter  results  for  selected 
3-dimensional  structures.  The  results  are  intended  to  verify  the  applicability  and  accuracy  of  the 
analysis  and  are  not  a  detailed  study  of  any  particular  structure. 


2  TECHNIQUE  DEVELOPMENT 


2.1  General  Procedure 

Consider  a  waveguide  circuit  element  made  up  of  a  number  of  cascaded  sections  as  shown  in 
Figure  2.1.  Each  section  is  uniform  in  the  direction  of  propagation  (chosen  to  be  the  z-direction), 
and  consists  of  a  number  of  connected  rectangular  regions.  Each  region  may  contain  an  arbitrary 
number  of  dielectric  layers.  A  brief  outline  of  a  transverse  resonance  mode  matching  technique  to 
determine  the  scattering  parameters  of  the  circuit  element  follows: 
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( i )  Subdivision  of  waveguide  element 

Each  section  (uniform  in  the  direction  of  propagation)  in  the  element  is  sub-divided  into  layers 
and  regions  comprising  homogeneous  and  isotropic  segments  of  rectangular  cross-section.  Layers 
aligned  vertically  and  not  separated  by  conducting  surfaces  are  combined  vertically  to  form 
subregions.  Vertically  aligned  subregions  are  in  turn  combined  to  form  regions.  Regions  are 
stacked  horizontally  to  construct  sections.  For  much  of  the  discussion  of  this  chapter,  only  simple 
regions  such  as  those  shown  in  Figure  2.1  which  consist  of  a  single  subregion  need  be  considered. 
(An  example  of  a  subregion  is  shown  in  Figure  2.5.)  Once  the  waveguide  element  is  subdivided, 
the  analysis  proceeds  by  first  analysing  all  of  the  subregions  of  a  section. 


( i  i )  Analysis  of  section  subregions 

The  transverse  resonance  technique  is  first  applied  in  the  y-direction  to  obtain  mode  eigenvalues 
and  eigenfunctions  for  each  subregion  in  each  section.  The  subregion  eigenmodes  propagate  in  the 
u-direction,  which  lies  in  the  x-z  plane.  This  concept  is  described  further  in  [12].  These 
eigenmodes  are  therefore  referred  to  as  w-modes.  For  a  given  free-space  wavenumber  kg,  the 
propagation  constant  k^,  or  the  w-mode  eigenvalue,  is  given  by: 


Note  that  while  ky  is  dependent  on  each  layer's  dielectric  constant  Ef  and  the  layer's  thickness, 
kx  (and  hence  k^)  is  constant  across  a  subregion  because  of  field  continuity  at  layer  interfaces. 
Field  continuity  at  subregion  interfaces  requires  kz  to  be  the  uniform  across  the  entire  section. 
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Therefore,  by  determining  the  subregion  eigenmodes  and  eigenfunctions,  the  effect  of  individual 
layers  is  determined,  and  the  analysis  proceeds  by  combining  subregions. 

(iii)  Analysis  of  individual  sections 

To  obtain  the  required  z-mode  eigenvalues  and  eigenfunctions  for  each  section,  the  transverse 
resonance  technique  is  applied  in  the  ;c-direction.  This  stage  of  the  calculation  is  generally  the 
most  computationally  demanding,  since  there  are  a  number  of  steps  involved  in  this  process 
which  must  be  repeated  many  times.  These  steps  are: 

(a)  Assume  an  initial  value  of  for  a  particular  z-mode  in  the  section. 

(b)  From  the  subregion  w-mode  data,  calculate  the  corresponding  x-mode  eigenvalues,  or 
propagation  constants  and  eigenfunctions  for  each  region.  The  propagation  constants 
are  given  by 

(c)  Apply  the  transverse  resonance  technique  in  the  a:-direction,  including  x-mode 
coupling  at  region  interfaces.  The  assumed  value  of  kz  is  varied  iteratively  and  steps  (b) 
and  (c)  repeated  until  the  transverse  resonance  condition  is  satisfied.  Such  values  of  k^ 
correspond  to  z-mode  eigenvalues  for  the  section.  The  z-mode  eigenfunctions  are  calculated 
from  x-mode  data. 

(iv)  Analysis  of  complete  element 

Having  obtained  the  field  expansion  in  the  individual  sections  expressed  by  the  z-mode 
eigenvalues  and  eigenfunctions,  the  scattering  parameters  of  both  the  sections  themselves  and 
the  discontinuities  arising  at  section  interfaces  can  be  determined.  By  cascading  these  scattering 
matrices,  the  scattering  matrix  for  the  entire  waveguide  element  is  readily  determined.  This 
constitutes  the  desired  output  from  the  analysis. 

2.2  Definitions  and  conventions 

The  rectangular  coordinate  system  from  [17]  shown  in  Figure  2.1  was  adopted,  with  unit  vectors 
denoted  by  (ux,  Uy,  Uz).  The  dependence  on  time  t  is  assumed  throughout  to  be  c/®*  where  j  -  V-1 
and  (0=27^  and  /is  the  operating  frequency. 

Let  the  total  electric  and  magnetic  field  in  the  waveguide  in  the  transverse  cross-section  of  the 
jth  section  be  given  by  vectors  E;  and  H/  respectively.  The  total  field  in  the  two-dimensional 
section  is  uniquely  described  by  the  tangential  fields  alone  so  we  may  let 

Ei=E,^+E,^.+E„ 

=  E,,+E,, 

where  Et;-  is  the  tangential  electric  field  and  similarly  for  the  magnetic  field.  If  the  field  is 
represented  by  a  sum  of  normal  modes  in  each  subregion  of  the  cross-section,  then  the  total 
tangential  fields  may  be  expressed  by  a  weighted  sum  over  all  possible  modes 
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i.e.  for  each  mode 


E„  =  S  vfer 


m=0 


H'V  iffiuin 

t.  =  2. 


m-Q 


E?  =  vrer 


where  and  are  equivalent  voltages  and  currents  which  are  functions  of  the  propagation 

coordinate  variable  only  and  e^  and  are  mode  functions  which  depend  on  the  transverse 
variables  (and  physical  geometry)  only.  These  equivalent  voltages  and  currents  can  be  used  in  a 
circuit  model  of  the  waveguide  e.g.  the  modal  transfer  admittance  between  mode  m  at  port  i  and 
mode  p  at  port  j  is  given  by 


tnp  ^h_ 

vP 

J 


2.3  Method 

Following  the  general  procedure  outlined  in  Section  2.1,  the  method  is  developed  in  greater 
detail.  The  transverse  resonance  technique  is  applied  first  in  the  y-direction  in  the  subregion 
analysis  and  then  again  in  the  3:-direction  in  the  section  analysis.  Once  the  required  number  of 
section  modes  have  been  calculated  for  each  section  using  the  transverse  resonance  approach,  a 
mode  matching  technique  is  used  to  calculate  the  scattering  parameters  of  the  complete  element 
formed  by  cascading  its  constituent  sections. 

2.3.1  Multilayer  parallel  plate  analysis  (subregions) 

The  first  step  is  to  apply  the  transverse  resonance  technique  in  the  y-direction  to 
determine  the  subregion  u-mode  eigenvalues  and  eigenfunctions.  This  mode  data  will 
implicitly  contain  the  effects  associated  with  individual  layers  in  a  given  subregion, 
reducing  the  required  analysis  to  the  subregion  level  for  each  two-dimensional  section. 

Consider  propagation  in  a  multi-dielectric-layer  parallel-plate  transmission-line 
structure  shown  in  Figure  2.2  [cf  [12],  p.  847].  This  is  a  classical  Sturm-Liouville  eigenvalue 
problem  [18].  To  find  the  eigenvalues  (propagation  constants)  and  eigenfunctions  (mode 
functions)  for  propagation  perpendicular  to  the  plates  (y-direction),  we  apply  the 
transverse  resonance  technique  in  this  direction.  Each  layer  can  be  then  considered  as  a 
portion  of  a  lossless  transmission  line  with  modes  propagating  in  the  y  direction,  as  shown 
in  Figure  2.2.  These  modes  are  referred  to  as  y-modes  because  of  their  direction  of 
propagation.  The  voltages  and  currents  on  the  transmission  lines  are  directly  related  to 
the  mode  functions,  and  hence  the  y-mode  electromagnetic  fields. 

Normally,  the  choice  of  modes  would  be  transverse  electric  (TE)  and  transverse  magnetic 
(TM)  to  the  direction  of  propagation  (y  in  this  case).  However,  alternative  representations 
using  modes  that  are  either  TE-to-x  and  TM-to-x  or  TE-to-z  and  TM-to-zmay  lead  to 
simpler  (and  computationally  more  efficient)  representations  under  particular  conditions 
such  as  subregions  or  sections  with  uniform  dielectric  constants.  The  choice  of 
representation  is  discussed  further  in  Section  2.3.2. 
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Figure  2.2  Equivalent  transmission  line  network  for  y-mode  analysis 

TE  and  TM-to-y  modes  are  used  in  multilayered  dielectric  subregions  since  there  is  no 
coupling  between  TE  and  TM  modes  at  dielectric  discontinuities,  resulting  in  a  more 
straightforward  analysis.  For  this  choice  of  eigenmode,  the  characteristic  admittance  of 
layer  I  with  relative  dielectric  constant  e;  is  given  by 


‘  mo 

y"  _  (OEqEi 

s 

where  (O  =  27tf  where  f  is  the  frequency,  po  ^0  are  the  permittivity  and  permeability 
respectively  of  a  vacuum. 

By  contrast,  TE  and  TM-to-x  and  z  modes  are  used  only  in  uniform  subregions,  which  consist 
of  a  single  dielectric  layer,  since  TE-TM  mode  coupling  would  otherwise  occur  at  dielectric 
discontinuities  between  layers.  For  imiform  subregions,  similar  admittances  can  be  defined 
for  TE,  TM-to-x  and  TE,  TM-to-z  modes  as  follows 

for  TE  -  to  -  X,  TE  -  to  -  z  modes. 


for  TM  -  to  -  X,  TM  -  to  -  z  modes 


v"  _  OJEq 

'~~r 

yi 


’  _  yi 


Y,= 


mo 


for  TE  -  to  -  y  modes 
for  TM  -  to  -  y  modes 
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The  restriction  that  these  modes  can  only  be  used  for  uniform  subregions  is  reflected  by  the 
fact  that  the  dielectric  constants  for  individual  layers  are  no  longer  contained  in  the 
modal  admittances. 

As  stated  previously,  the  propagation  constant  in  the  y-direction  for  the  layer,  is 
related  to  the  propagation  constant  ky  by 

^>7  ~  ^0^1  ~  ' 

The  subregion  analysis  proceeds  by  first  assuming  the  maximum  possible  value  of  ky.  The 
admittances  of  the  transmission  lines  shown  in  Figure  2.2  are  evaluated  for  each  y-mode. 
Using  these  admittance  values,  the  modal  voltages  and  currents  used  to  construct  the  mode 
eigenfunctions  are  calculated  at  the  interfaces  of  each  layer.  Depending  on  the  boundary 
condition  (electric  or  magnetic  wall),  either  the  y-mode  voltage  or  the  current  is  set  to  zero 
at  the  lower  plate.  If  the  value  of  ky  is  such  that  the  voltage  and  current  at  the  top  plate 
satisfy  the  required  boundary  condition,  then  ky  is  a  valid  wavenumber,  otherwise  the 
value  of  ky  is  iteratively  adjusted  until  the  boundary  condition  is  satisfied.  Note  that 
this  method  is  preferable  to  using  the  usual  method  of  finding  a  resonance  [4]  in  the  total 
admittance  found  by  summing  the  transmission  line  admittances  in  two  directions  at  an 
arbitrary  plane,  because  the  voltage  and  current  remain  finite  whereas  the  admittance 
may  approach  ±o°.  The  search  for  valid  mode  solutions  is  considerably  more  difficult  with 
the  presence  of  poles. 

The  y-mode  voltage  and  current  at  the  end  of  the  line  representing  the  layer  are 
obtained  from  the  voltage  and  current  at  the  start  of  the  line  by 

Vi+i  =  coskyid  Vi  -  —sinkyid  Ii 

Il+l  =  coskyid  I  I  -  jYi  sinkyid  V/. 

with  the  signs  consistent  with  the  convention  shown  in  Figure  2.2. 

For  a  given  y-mode,  the  admittance  matrix  for  a  length,  d,  of  transmission  line 
representing  layer  I  is 

-7/  (coth  ayid+ j  cot 7/  (cosech  Uyid+j  cosec  Pyid^ 

7/(cosech  ttyid+j  cosec)3yd)-7;(coth  a yid+ j  cot Pyid^ 
where  Y/  is  the  modal  admittance  of  the  layer  for  kyi  -  |3yj  +  jOy;. 

Using  the  appropriate  admittance,  the  transmission  line  parameters  ^y{y)  Iy(y)  are 
used  to  construct  mode  functions  which  represent  the  electric  and  magnetic  fields  of 
TE-to-y,  TM-to-x,  and  TM-to-z  modes.  Similarly,  the  parameters  Vy(y)  and  Iy(y)  are  used 
to  construct  mode  functions  which  represent  the  electric  and  magnetic  fields  of  the 
TM-to-y,  TE-to-x,  and  TE-to-z  modes.  Because  these  parameters  are  used  to  represent 
propagation  in  the  y-direction,  they  are  termed  the  y-mode  voltages  and  currents. 

The  amplitudes  of  the  transmission  line  voltages  and  currents  used  to  define  the  y-mode 
mode  fimctions  are  normalised  to  yield  the  following  orthonormality  relations: 

\'^'y"y{y)yy(.y)*dy  =  d,j 
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for  TE-to-y,  TM-to-x,  and  TM-to-z  modes. 


p+i 

■'yi 


Iyiy)iy(yy 

£(y) 


dy  =  5y 


for  TM-to-y  modes  and 


i'yi.y)iy{yVdy  =  5.. 


for  TE-to-x  and  TE-to-z  modes.  In  the  above  expressions,  8m-  is  the  Kronecker  delta. 

If  •' 

denotes  the  complex  conjugate,  and  y/  and  yi+i  are  the  lower  and  upper  edges  of  the 
layer.  The  phase  ambiguity  in  the  normalisation's  allows  the  choice  of  V  (y)  to  be 

y 

always  real,  and  I  (y)  to  be  always  imaginary. 

y 


The  orthonormality  of  Vy(y)  and  Iy(y)  are  used  in  Section  13.1.1  to  calculate  the  mode 
coupling  between  adjacent  regions  in  order  to  represent  the  entire  two-dimensional  cross 
section. 


2.3.2  Two  Dimensional  (Section)  Analysis 

Section  2.3.1  has  described  the  subregion  analysis  which  is  used  to  reduce  the  problem  of 
analysing  the  various  x-y  cross-sections  of  the  structure  to  the  subregion  level.  The  effects 
of  individual  layers  are  implictly  included  in  the  y-mode  and  w-mode  propagation 
constants  and  mode  functions.  The  remainder  of  the  two-dimensional  section  analysis  is 
concerned  with  combining  subregions  into  sections. 

The  two-dimensional  analysis  described  in  this  section  can  be  applied  to  a  wide  variety  of 
waveguiding  structures  including  dielectric  guides,  finline,  microstrip,  and  coplanar 
waveguide.  The  results  of  the  analysis  of  a  range  of  waveguides  are  presented  in  Section 
3.1.  The  z-mode  propagation  constants,  kz,  can  be  determined,  along  with  the  modal 
fields.  These  fields  can  be  used  to  calculate  user-defined  z-mode  impedances,  and  to 
calculate  the  z-mode  coupling  between  connected  sections.  This  coupling  is  used  to  cascade 
the  sections  to  construct  the  entire  waveguide  element. 

To  fully  describe  the  fields  for  modes  propagating  in  the  z-direction  in  a  given  section 
(z-modes),  the  transverse  resonance  technique  is  applied  to  an  equivalent  transmission  line 
network  representing  propagation  in  the  x  direction.  The  modes  represented  by  this 
network  are  referred  to  as  x-modes  because  of  their  propagation  direction.  The  previously 
determined  y-mode  fields  which  were  determined  using  transverse  resonance  in  the  y- 
direction  are  combined  as  weighted  sums  to  generate  the  fields  of  each  z-mode  in  the 

section.  The  propagation  constants  of  the  z-modes  are  {m  =  0,  1,  ...,  Nz),  where  N2  is 
the  total  number  of  z-modes  to  be  considered.  Ideally,  an  infinite  number  of  modes  exist, 
but  this  number  must  be  truncated  for  practical  purposes. 

The  x-mode  voltages  and  currents  in  the  x-mode  network  defining  the  contribution  of  the 
y-mode  fields  to  the  total  fields  of  the  z-mode  are  and  1^”.  Figure  2.3  shows 
the  x-mode  transmission  line  network  for  a  simple  ridged  waveguide  cross-section.  Note 
that  in  general,  the  transmission  lines  representing  individual  x-modes  in  each  region  are 
now  coupled  at  region  interfaces.  The  boundary  conditions  at  the  sides  of  the  cross-section 
determine  how  the  x-mode  transmission  lines  are  terminated.  The  terminations  are  short 
circuits  for  an  electric  wall,  open  circuits  for  a  magnetic  wall,  and  the  characteristic  modal 
admittances  (match)  for  an  open  boundary  such  as  in  a  parallel-plate  transmission  line. 
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Figure  2.3  Ridged  waveguide  cross-section  showing  ar-mode  transmission  line  network 


The  X-  mode  and  y-mode  voltages  and  currents  are  combined  as  follows  to  represent  the 
components  of  the  fields  transverse  to  the  x-direction 

E  =  E^Ux  +  EyViy  +  E^n^ 

=  Ex^ix  + 

where 

^r  =  SS''rWe?-(>’.z) 

m  n 

and 

ll  =  Hx\ix  +  Hy\iy^H^n^ 

where 

m  n 

Ny  and  Nz  are  the  number  of  transverse  eigenmodes  and  the  number  of  z-modes 
respectively.  The  modal  fields  and  hj  can  be  derived  in  terms  of  the  y-mode  voltages 
and  currents  using  the  expressions  derived  in  Appendix  I.  The  three  choices  of  eigenmodes 
discussed  in  Section  2.3.1,  ie,  TE  and  TM  to  either  x,  y,  or  z,  have  features  relevant  to  a 
discussion  of  the  section  analysis.  The  facility  to  choose  from  among  the  three  possible 
types  of  mode  functions  is  useful  to  characterise  the  z-modes  in  the  structure,  to  allow 
various  types  of  excitation  to  be  applied  in  a  straightforward  fashion,  and  to  improve 
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efficiency.  In  particular,  each  choice  of  eigenmode  type  has  specific  advantages  in  special 
circumstances.  These  advantages  include: 

(i)  For  regions  with  layers  of  different  dielectric  constants  (i.e.  dielectric  interfaces  in 
the  horizontal  or  x-z  plane),  TE-to-y  and  TM-to-y  modes  must  be  used.  Of  the  three 
choices,  only  these  modes  have  the  property  that  no  coupling  occurs  between  different 
modes  at  dielectric  interfaces,  simplifying  the  subregion  analysis  of  Section  2.3.1  for 
multilayered  subregions.  The  lack  of  coupling  arises  from  the  orthogonality  condition  in 
the  definition  of  the  modes. 

(ii)  To  simplify  the  analysis  of  sections  containing  dielectric  slabs  with  interfaces  in  the 
vertical  or  y-z  plane,  TE-to-x  and  TM-to-o:  modes  can  be  used,  since  no  coupling  occurs 
between  the  TE  and  TM  modes  at  dielectric  interfaces  that  are  normal  to  the  reference 
direction  for  the  modes.  As  discussed  in  Section  2.3.1,  the  subregions  in  the  section  must 
uniform,  otherwise  there  will  be  coupling  between  TE  and  TM-to-x  y-modes  at  layer 
interfaces.  Coupling  between  mode  functions  in  different  subregions  results  only  from 
discontinuities  in  the  conducting  boundary  of  the  section  rather  than  from  changes  in 
dielectric  constants,  making  these  modes  an  ideal  choice  for  selected  dielectric 
waveguides. 

(iii)  For  sections  filled  with  a  medium  of  uniform  dielectric  constant,  TE  and  TM-to-z 
modes  are  uncoupled.  Therefore,  this  choice  of  modes  can  lead  to  considerable  savings  in 
computational  overhead  arising  from  simplified  mode  coupling  between  subregions.  In 
addition,  the  z-modes  subsequently  calculated  using  the  x  and  y-modes  in  each  sections  can 
be  readily  identified  as  TE  or  TM-to-z,  which  can  simplify  the  physical  interpretation  of 
the  results. 

A  number  of  two-dimensional  transverse  resonance  analyses  have  been  developed  including 
those  described  in  [12],  [13],  and  [19]-[25].  Although  varying  significantly  in  detail,  the 
general  procedure  is  similar  in  all  cases. 

2. 3.2. 1  Mode  Function  Derivation 

In  this  section,  the  components  of  x-mode  fields  transverse  to  the  x-direction  are  related  to 
the  voltages  and  currents  on  the  y-mode  transmission  lines.  This  step  is  essential  for 
combining  subregions  into  a  complete  two-dimensional  section  using  the  transmission  line 
network  discussed  in  Section  2.3.2. 

The  general  technique  used  to  generate  the  mode  functions  follows  Felsen  and 
Marcuvitz  [10],  and  is  derived  in  Appendix  I.  The  derivation  of  TE  and  TM-to-y  mode 
functions  is  presented  in  this  section,  with  the  corresponding  derivations  for  TE,  TM-to-x 
and  TE,  TM-to-z  modes  included  in  Appendix  II  for  completeness. 

TE-to-y  modes 

Since  ey  =  0,  the  modal  fields  transverse  to  x  are 

ej{y.z)  =  e^(y)e~^'‘z'u^ 


where  s^(y),  h^{y)  and  h^{y)  are  field  coefficients  which  depend  on  y  only  and  the  prime 
indicates  a  TE  mode.  Following  the  derivation  in  Appendix  I,  the  transverse  components  of 
the  modal  electric  and  magnetic  fields  are  related  by 


k^Y^hj{y,z)  =  o)e 


1 


’UxXerM 
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Writing  the  field  components  explicitly  and  cancelling  the  e'i'^z^  factor  yields 


(oe[{  3^  ,  ,2 


Using 


9  ,  ,.2  1  '  .  .,2  ,2x  ' 


=  jtV 

and  choosing  the  arbitrary  admittance  factor  to  be 

F  = 

^  k^CO/l 

then  the  following  relationships  between  the  electric  and  magnetic  field  components 
transverse  to  x  emerge 

hy  —  —e^ 


jk.  de, 

kl  ay 

The  field  quantities  can  be  related  to  the  voltages  and  currents  on  an  equivalent 
transmission  line  network  as  shown  in  Figure  2.2  by  setting 

4  =  ^y(F) 

This  allows  the  magnetic  field  components  to  be  specified  as 


hy  =  -Vy(y) 


jk.  d  . 
tidy 

=  Kh.±Uy) 


Y  y 
K  ^0 


since,  according  to  transmission  line  theory. 


^  y'yiy)  =  -jky  4-  iy{y)  =  -jko  Iy(y) 
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TM-to-y 

Proceeding  in  a  manner  identical  to  the  TE-to-y  mode  treatment  and  using  hy  =  0,  the 
modal  fields  transverse  to  x  are 

hT()'.z)  = 

eT(y,z)  =  ey{y)e~^^z^ny  +  e^{y)e~^^z^v^^ 

where  h^(y),  e^{y)  and  e^{y)  are  field  coefficients  which  depend  on  y  only  and  the 
double  prime  indicates  a  TM  mode.  Following  the  derivation  in  Appendix  I,  the  transverse 
components  of  the  modal  electric  and  magnetic  fields  are  related  by 

kxZx^T{y>^)  =  ^  ‘hrMxUx 

L  <y^  £(y)  J 

Writing  the  field  components  explicitly  and  cancelling  the  e'j^z^  factor 


o)ti  ^  1  a  , .  2 


k,Z,ieyny+e,u^)-  ^  N  jK 


1  dhy 


"  ^  e/(y)ay  ^ 


A  simplification  can  be  made  using 


3  1  3  , 2  I  3  1  dK  . 2  K 

- — — +  ^n  - —^  +  kti—^ 


^y  £i(y)  ^y  ^  J  ay  e;(y)  ay  e/(y) 


■(^0^/(y)~^yz) 


’"leiV) 


since  e;(y)  is  constant  within  each  layer.  By  choosing  the  arbitrary  impedance  factor  to  be 


kx(0£Q 

then  the  following  relationships  between  the  transverse  electric  and  magnetic  field 
components  emerge 

(t 

y 


/=zAl  ±  ^ 

^  kl  e,  3y 

The  field  quantities  can  be  related  to  the  voltages  and  currents  on  the  equivalent 
transmission  line  network  by  setting 


K  =  iyiy) 


This  allows  the  electric  field  components  to  be  specified  as 
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and 


It 


Iy(y) 

£1 


_  -JK. 


■  -2Z 


u 

kJcQ 


\  d  " 
—Yh^y) 

El  dy 

■yQVyiy) 


2.3.2.2  Mode  Coupling 

The  x-mode  fields  inside  each  subregion  can  now  be  expressed  in  terms  of  y-mode 
transmission-line  voltages  and  currents.  To  complete  the  analysis  of  a  complete  section, 
subregions  are  stacked  together  horizontally  and  the  coupling  between  :r-modes  in  adjacent 
subregions  is  evaluated.  This  coupling  satisfies  the  field  continuity  across  the  subregion 
interface  shown  in  Figure  2.4.  The  evaluation  of  this  x-mode  coupling  at  subregion 
interfaces  is  discussed  in  this  section. 

Mode  functions  in  adjacent  subregions  will  be  coupled  at  discontinuities  in  the  conducting 
region  boundaries  or  dielectric  constants  between  the  two  subregions.  This  coupling  can  be 
evaluated  from  continuity  of  the  transverse  electric  and  magnetic  field  components  in  the 
y-z  plane  of  the  discontinuity,  and  from  the  boundary  conditions.  Consider  an  interface 
between  simple  regions  each  consisting  of  a  single  subregion,  where  a  change  in  the  region 
heights  or  dielectric  constants  occurs  as  shown  in  Figure  2.4. 
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From  the  definitions  of  the  x-mode  voltages  and  currents  stated  previously  in  the 
introduction  to  Section  2.3.2,  the  field  components  in  region  A  tangential  to  the 
discontinuity  are  given  by 

Er,  =SIV'";We7-/>’.z) 

m  n 

Hr, 

m  n 

and  similarly  for  region  B 

N^Nf 
m  n 

m  n 

where  the  summation  is  over  both  TE  and  TM  eigenmodes.  The  ratio  between  the  number  of 
y-modes  in  regions  A  and  B  is  determined  by  the  relative  convergence  criterion 

<  h 

Relative  convergence  and  its  explicit  role  in  the  solution  of  characteristic  equations 
encountered  in  general  mode  matching  analyses  is  discussed  in  detail  by  Leroy  [26]. 

The  field  matching  at  the  discontinuity  requires  that 

g  _  ,  over  aperture 

lo,  otherwise 


and 


over  aperture. 

By  equating  the  modal  expansions  for  the  fields  in  regions  A  and  B  in  the  aperture,  the 
coupling  relations  between  the  x-mode  voltages  V^^(x)  and  v'g^(x),  and  between  the  x- 

mode  currents  and  l^^(x)  can  be  evaluated  for  the  m*^  z  mode  for  all  m.  This 

procedure  was  followed  to  determine  the  coupling  relations  for  each  of  the  three  types  of 
basis  modes,  ie.,  TE,  TM  to  x,  y,  and  z.  The  coupling  expressions  were  adjusted  into  a  form 
that  guaranteed  conservation  of  complex  power  across  the  discontinuity  for  a  finite  number 
of  mode  functions  in  each  region,  while  still  maintaining  the  field  matching  for  the  ideal 
case  where  an  infinite  number  of  mode  functions  are  used.  In  the  following  section,  the  TE 
and  TM-to-y  mode  function  coupling  expressions  are  derived.  The  procedure  is  to  first 
match  the  two  transverse  components  of  the  electric  field  (the  y-component  and  the  z- 
component),  using  the  appropriate  orthonormality  relations  to  simplify  the  expressions 
and  then  repeat  the  procedure  for  the  two  transverse  components  of  the  magnetic  field.  The 
corresponding  derivation  for  TE,  TM-to-x  and  TE,  TM-to-z-modes  is  given  in  Appendbc  III. 
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Electric  Field  Matching 
y-component 

fjTM  ""  ,  .  f^TM  "P  ,  , 

B  mn  I  By  Y  ^ 

„  £Biy)  ~  p  ^A(y) 

„n  ^ 

After  multiplying  both  sides  by  Ig  (y)  and  integrating  along  the  region  interface 
extending  from  y;  to  y^  =  y;  +  IzA  shown  in  Figure  2.4,  the  orthonormality  relation 

allows  the  coupling  between  the  x-mode  TM  voltages  to  be  determined 


S...,'f4'WCW 


P  y, 


^A^y) 


where 


Ju 

A„„  = 

np  J 


e.(y) 


z-component 


N™ 

uTtl, 


^  <mn  <n  ^  Ir'’^  (n£[\  '""P  "P 

S  %  Vs<y)-  S  Fb,  Vs(y) 


’’  ^“pb 

M™ 


.mi  ,i  ^'A  •  • 

=  I  ^A,^A,(y)-  £  ^A/y) 


j 


After  multiplying  both  sides  by  Vg  (y)*  and  integrating  along  the  region  interface,  the 

y 

orthonormality  relation 

|Q^Vy)^is/y)*^y=^y 

allows  the  coupling  between  the  x-mode  TE  and  TM  voltages  to  be  determined 

N™  y 

,mn  B  u  *  "P 

^B,  -  S  J  ^Byiy)  yBy(y)dy 


.,2  '  B, 

P 

nI^  .yu 

^  .mi  r  .n  *,  •' 


k^COEn  *  "j 

V  _0  T/  f  T/_  r^.v  T/  . 


=  S  J  ysMyliy^dy-  vliyfv^y^y 

i  yi  j  ^'^ja  yi 
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which  can  be  written  using  suitably  defined  matrices  F,  C,  and  D  as 


where 


Jc  —  ,,2 

np  .  "2 


“pb  yi 


C.i-lv^WK/rtrfy 


k^coe,, 


D^-^pT^\vl(yrviydy 


"Ja  yi 


In  anticipation  of  later  modifications  to  ensure  complex  power  conservation  across  the 
region  interface,  a  matrix  O  is  defined  by 


v=^jv;wv;w4y 


’^pa  yi 


Magnetic  Field  Matching 
y-component 


Nf 


f  mn  I  n  *  mp  \  p 

£  I  A,  VA,(y)=  £  VB,(-y) 


nV 


mp  I  p 

I 

X 


,11  ^ 

After  multiplying  both  sides  by  (y)  and  integrating  along  the  region  interface,  the 
orthonormality  relation 


Jy,  y  y  A 


allows  the  coupling  between  the  x-mode  TE  currents  to  be  determined 

N?  yu 


I  mn  «  1  tmp  r  <11  *  >  P 

h.  =£v  v'./yW)’ 


P  yi 


=£c.-/,: 


where 


=  jv'lSyKfy^dy 
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z-component 


iim«  «n  .Ay  k'^WLln  ”^P  'P 

I  h,  U)l^^(y)+  S  MlA,(y) 


p  I^Up^ 

7,/n, 


■A  "fni  kiPwLln  '"‘i  '■' 

=  S  /b,  (A:)/B,<y)+  I  (x)IB: 


i 


After  multiplying  both  sides  by  (y)*  and  integrating  along  the  region  interface,  the 

y 

orthonormality  relation 

ch  1  Jk  ,, 

[  (y) - (y)  = 

Jo  £j^(y)  y 

allows  the  coupling  between  the  x-mode  TE  and  TM  currents  to  be  determined  as 

Y£  Mrt  *  '  P 

’a,  +  L  X  J 

p  *“ba  ^ 


^  ^  £A(y)  ^ir  ^  J  e^Cy) 


j  yt 


which  can  be  written  using  suitably  defined  matrices  F,  A  and  j8  i 


..  N™ 


,.mn  .A  <mp  ®  ^  i.mi  o  ,'"y 

I  A.  +  I  r„/A,  =  S  +  X  a/b, 

n  1  I 


where 


C  ! 

"■  C  I 

In  anticipation  of  later  modifications  to  ensure  complex  power  conservation  across  the 
region  interface,  define  a  matrix  G  is  defined  by 
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2.3.2.3  Complex  Power  Conservation 

The  condition  that  the  transverse  electric  and  magnetic  fields  be  exactly  matched  across 
region  interface  discontinuities  will  only  be  satisfied  in  the  ideal  case  where  an  infinite 
number  of  x-modes  are  used.  For  practical  purposes,  the  number  of  x-modes  must  be 
truncated  to  a  finite  number,  so  that  exact  field  matching  cannot  be  achieved  in  general. 
However,  a  physical  constraint  that  power  flowing  through  the  region  interfaces  must  be 
conserved  can  be  enforced  on  the  field  matching  for  an  arbitrary  finite  number  of  x-modes. 
That  is,  the  power  transmitted  in  the  x-direction  from  region  A  must  equal  the  power 
received  in  region  B,  and  vice  versa.  Both  the  real  and  imaginary  components  of  the  power 
flow  are  matched  at  the  interfaces  to  enforce  this  most  fundamental  physical  constraint  on 
the  x-mode  coupling  at  region  interfaces.  This  general  procedure  has  previously  been 
applied  to  a  two-dimensional  mode  matching  analysis  by  Mansour  and  MacPhie  [27], 
although  only  in  homogeneous  structures.  Omar  and  Schunemann  [28]  have  shown  that 
conservation  of  complex  power  can  be  applied  to  inhomogeneous  structures,  and  that  it 
follows  directly  from  field  matching.  This  conclusion  has  been  verified  by  this  work. 

By  constructing  column  vectors  for  the  x-mode  voltages  and  currents  for  the  m*^^  z-mode 
(with  the  TE  and  TM  components  separated),  and  using  the  matrices  defined  in 
Section  2.3.2.2,  the  coupling  expressions  from  the  field  matching  can  be  written  in  matrix 
form  as 


The  power  flow  through  the  region  interface  on  either  side  of  the  discontinuity  for  the  m*^ 
z-mode  is  determined  from  the  fields  in  region  A  and  region  B  using  the  x-mode  voltages 
and  currents  as  follows 
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The  coupling  between  the  modal  fields  in  regions  A  and  B  is  required  to  be  such  that  the 
power  coupled  from  region  A  to  region  B  is  the  same  as  the  power  coupled  from  region  B  to 
region  A.  For  region  A 


and  for  region  B 


So  to  ensure  conservation  of  complex  power  for  an  arbitrary  number  of  x-modes,  the 
following  substitutions  must  be  made 

P  -^G+A 

D  ->CO 

For  these  substitutions  to  have  physical  relevance,  they  must  be  automatically  satisfied 
for  an  infinite  number  of  modes,  showing  that  the  fields  on  each  side  of  the  discontinuity 
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will  still  match  if  infinite  modes  are  used.  This  can  be  shown  using  the  completeness  of  an 
infinite  set  of  eigenmodes  as  follows 


[GA]„p  =  S 


ttg  tfi  Jc  »»p 

kTcojdo  hy  ^y^hy  ^y^  .  ^y  ^ 


.2 


yi 


-dy  j 


yi 


£A(y') 


dy 


tfl  ♦  "P  TAA 

k^ma  'b.  <y)  'a.  (>■ )  h.  (’'■ ) 


z=-L 


a 


^Aiy') 


yi  yi  “  ■  ^ 

<  i 

jSM  as  required. 
mp 


£B(y) 


dydy' 


Similarly 


TE 

'«  *  '9  k^OJEn  '9  *  "P 

[c*i«,,=  E  J  ''b,w  I  %()'■) 

>-/ 


^  yi 


a-F 


1  ^;w’v;(y)S  ^I(y)*Cw<^</y 


“/^A  >/  y/ 


K,  COEr)  f  •«  *  "P 

=  -^T2-^J^fi  (>-) 

=  Dfjp  as  required. 


The  final  form  for  the  coupling  which  guarantees  conservation  of  complex  power  for  any 
number  of  x-modes  in  each  region  is  therefore 


2.3.2.4  Constructing  Admittance  Matrices 

So  far,  only  simple  region  interfaces  have  been  considered,  i.e.  regions  that  do  not  contain 
subregions  separated  by  conducting  boundaries.  At  this  point,  it  is  necessary  to  consider 
more  complex  regions  consisting  of  more  than  one  subregion  such  as  shown  in  Figure  2.5.  The 
procedure  is  to  calculate  the  two-port  x-mode  admittance  matrices  of  the  subregions  in  a 
given  region,  calculate  the  coupling  matrices  for  the  interfaces  between  connected 
subregions,  and  then  combine  the  admittance  matrices  and  the  coupling  to  produce  a  new 
admittance  matrix  for  the  region  which  now  includes  the  coupling  to  other  regions.  In  this 
manner,  a  series  of  multiport  admittance  matrices  can  be  constructed  to  represent  each 
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region  in  a  given  section.  When  cascaded  together  in  the  correct  manner,  these  admittances 
represent  the  entire  section  under  consideration.  Note  that  the  independent  subregions  in  a 
given  region  are  vertically  aligned  but  separated  by  conducting  metal  surfaces,  so  that  no 
y-mode  coupling  exists  between  them. 


region  A  region  B 


Figure  2.5  Complex  region  interface  with  independent  layers 

The  coupling  expressions  relating  the  x-mode  voltages  and  currents  for  the  m*^^  z-mode 
between  connected  subregions  as  derived  in  Section  2.3.2.3  are  used  in  the  form 


for  "step-up"  transitions  (height  of  subregion  A  smaller  than  height  of  subregion  B),  and 

0  j  Pv  l"!^/ 

.«J '  L  p.  !  0  Jw.- 

for  "step-down"  transitions  (height  of  subregion  A  larger  than  height  of  subregion  B).  Pi 
and  Py  are  the  current  and  voltage  couplings,  respectively.  These  coupling  expressions  are 
listed  in  Appendix  IV  for  the  various  choices  of  eigenmode  types  in  subregions  A  and  B. 

The  coupling  expressions  above  show  that  currents  in  the  smaller  subregion  can  be 
expressed  in  terms  of  currents  in  the  larger  subregion,  but  not  vice  versa.  This  is  because  of 
the  lack  of  a  suitable  constraint  on  the  tangential  component  of  the  magnetic  field  of  the 
larger  subregion  on  the  aperture  wall.  Since  magnetic  field  matching  can  only  be 
considered  over  the  aperture  itself,  it  is  only  possible  to  determine  the  currents  in  the 
smaller  subregion  from  the  currents  in  the  larger  region.  To  determine  the  currents  in  the 
larger  region  from  the  currents  in  the  smaller  region  requires  the  nonexistent  magnetic  field 
boundary  condition  on  the  aperture  wall. 

Such  a  boundary  condition  does  exist  for  the  electric  field  however,  namely  that  the 
tangential  electric  field  is  required  to  be  zero  on  the  aperture  wall.  This  condition  is  used 
to  derive  the  voltage  coupling  relating  the  voltages  in  the  larger  subregion  to  the  voltages 
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in  the  smaller  subregion.  However,  if  the  tangential  electric  field  must  be  zero  on  the 
aperture  wall  then  an  arbitrary  electric  field  is  not  allowed  in  the  larger  subregion. 
Therefore,  it  follows  that  a  coupling  expression  relating  the  voltage  in  the  smaller 
subregion  to  an  arbitrary  voltage  in  the  larger  subregion  cannot  be  defined.  These 
considerations  are  also  discussed  in  [28].  Note  that  it  is  possible  to  numerically  invert  the 
coupling  matrices  to  circumvent  this  problem  if  the  same  number  of  x-modes  are  used  in 
each  subregion.  However,  not  only  does  this  approach  violate  relative  convergence,  it  also 
leads  to  a  very  ill-conditioned  formulation. 

A  subregion  with  a  two-port  modal  admittance  matrix  Ys  contained  in  a  given  region 
which  is  coupled  to  M  subregions  in  a  region  on  the  input  (left)  side  and  N  subregions  in 
another  region  on  the  output  (right)  side  is  shown  schematically  in  Figure  2.6.  The 

voltage  and  current  couplings  at  the  m*^  input  port  are  written  as  Py  and 

respectively,  and  the  corresponding  couplings  at  the  n^^  output  port  as  Py^  and  F”^ 
respectively. 


Figure  2.6  Multiport  admittance  block  constructed  from  two-port  subregion  coupled  to  other  subregions 

The  voltages  at  the  parallel  outputs  of  the  coupling  ports  sum  to  yield  the  total  voltage  at 
the  subregion  admittance  ports.  Therefore,  using  the  x-mode  voltage  and  current  coupling 
expressions,  the  voltages  and  currents  at  the  input  and  output  ports  of  the  subregion  were 
given  by 

M 

=^Vy  for "  step  -  up"  transition,  undefined  otherwise. 

m 
N 

=  ^Py  V”  for "  step  -  down"  transition,  undefined  otherwise. 

n 

jm  _  for "  step  -  down"  transition,  undefined  otherwise. 

/”  =  /2  for  "step  -  up"  transition,  undefined  otherwise. 

The  restrictions  discussed  previously  relating  to  the  boundary  transitions  that  may  be 
represented  by  the  coupling  expressions  do  not,  in  practice,  restrict  the  geometry  of 


24 


DSTO-RR-0027 


structures  that  can  be  analysed.  Any  arbitrary  cross-section  can  be  divided  such  that  all 
the  region  couplings  are  included  at  the  input  or  output  ports  of  the  subregion  admittances 
as  required.  Step-down  transitions  are  assigned  to  the  output  of  the  adjacent  subregion  to 
the  left,  while  step-up  transitions  are  assigned  to  the  input  of  the  adjacent  subregion  to  the 
right. 

The  current  at  the  input  and  output  ports  of  the  subregion,  Ig  and  Ig,  can  be  expressed  in 

terms  of  the  input  and  output  voltages,  Vg  and  Vg  using  the  two-port  subregion  admittance 
as 


4  =  Yo  V‘s  +  Y.  VI 
a  Ooo  oqI 


Therefore,  the  current  at  the  m*^  input  port  can  be  written  as 

M  N 

n  p 

Similarly,  the  current  at  the  n*^^  output  port  is  given  by 

M  N 

C  =  XPvVf + PJAs,, 

n  p 

Therefore,  a  multiport  admittance  matrix  which  represents  both  the  subregion  and  the 
coupling  to  regions  on  each  side  containing  multiple  subregions  can  be  defined 


-vr 

[YmoqJ  I 

v'm 

1— 1  j 
O  o  1 

v° 

Y  O 

L  tN 

—  V  — 

^  0 

where 


''MS  =  P?,Ys„Pf, 

YMS?  =  P?T5i,Pf^ 

YMrg=P?„Ys,„Pf, 

YM?f=P?Y.  Pf_ 


2.3.2.5  Cascading  Admittance  Matrices 

The  admittance  matrices  derived  in  Section  2.3.2.4  which  contain  both  the  subregion 
modal  admittances  and  the  associated  couplings  are  cascaded  to  represent  an  entire 
section.  The  cascading  of  two  multiple-port  admittances  Ya  and  Yg  to  form  a  single 
multiple-port  admittance,  Yc,  as  shown  in  Figure  2.7,  is  similar  to  simple  2-port 
admittance  cascading. 
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Figure  2.7  Multiport  admittance  cascade 

The  currents  at  the  input  and  output  ports  for  admittance  Y\  are  related  to  the  voltages  by 

N  N' 

4=IvaS^+IYa5?v, 

P  9 

N  N' 

p  q 

and  similarly  for  Yb 

-4=iYB;sv'p+iYBSv; 

p  q 

tr  N" 

4=EYBf?Y;,+XYBj?Y; 

p  9 

In  matrix  form 


i  =  YaooV+YaoiV 

r  =  YaioV + YauV  =  -[yboqV’  +YboiV"' 

I"=YbioV+YbiiV" 

Substituting  for  I' 

Ya„+Yboo]v  =  -[Ya,oV  +  Ybo,V"] 

Defining 

Yinv  “  [Yaii  Yroo 
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and  substituting  for  V  allows  I,  and  I"  to  be  determined  in  terms  of  V,  and  V" 

I  =  [Yaoo  -  YAoiYi„vYAio]Y -  YAoiYinvYBoiV" 

r= -YfiioYinvYAioV  +  [Ybh  -  YBioYinvYBoijv" 

In  matrix  form 

YCOO  I  ''^COI 


''^Cio  I  Yen. 

[Yaoo  "  YaoI  Yinv  YaiqI  |  "  Y^gi  Yinv  Ybqj 


■  Ybio  Yinv  Yaio  |  IYbh  -  Ybiq  Yinv  YboiI 

2.3.2.6  Junction  Admittance  Calculation 

By  applying  the  multiport  admittance  cascading  discussed  in  Section  2.3.2.5  in  an 
appropriate  manner,  the  section  to  be  analysed  can  be  reduced  to  two  multiport 
admittances  Yl  and  Yr  representing  the  structure  on  the  left  and  right  sides  of  a  junction 
plane.  The  position  of  the  junction  plane  can  be  arbitrary,  but  is  usually  set  at  the  left  side 
of  the  region  containing  the  subregion  with  the  minimum  number  of  y-modes.  This 
subregion  is  termed  the  junction  subregion.  Selecting  the  subregion  with  the  smallest 
number  of  y-modes  improves  computational  efficiency.  This  situation  is  shown 
diagrammatically  in  Figure  2.8. 

A  two-port  admittance  matrix  Yj  for  the  junction  was  determined  in  order  to  apply 
transverse  resonance  to  the  two-dimensional  section.  This  two-port  is  constructed  from  the 
admittances  "seen"  by  looking  to  the  left  and  right  of  the  junction  plane.  By  definition,  the 
junction  subregion  is  at  port  nj  of  Yr  and  port  nj  of  Y^.  From  the  definition  of  Y-parameters, 
the  current  at  the  n^  port  of  Yl  is  given  by 

P 

and  similarly  the  current  at  the  n^^  port  of  Yr  is 


N 


P 


Connecting  all  the  ports  of  Yj  and  Yr  except  for  the  jimction  port,  nj,  leads  to 

/2  =  -/"  and  forn^tnj. 
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Figure  2.8  Reduction  of  cross  section  to  2-port  junction  admittance 

The  input  and  output  currents  for  the  junction  port,  and  ,  are  given  in  terms  of  the 
voltages  at  the  other  ports  by 

p 

P*nj 

p^nj 


p^nj 
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Now  consider  the  currents  at  the  n*^  port  of  Yi  for  n7!:nj 

ptnj 


and  sinailarly  for  Yr 


N 


Y,y7vI 

P*nj 


Therefore,  since  /£  —  ~Ir  for  n?snj. 


N 


liyT+yfw^R 

p^nj 


Define 


where 


Yi„v=A 


-1 


Kp=\yT 


so  that  for  n^^nj 


N 


p^nj 


Hence,  substituting  for  Vr  (n?itnj) 


'l‘  = 


N  N 


vl'- 


i  lyl'^yS^jr 


L 

p:^nj  qi^nj 

J 

lP*njq:>tnj  J 

N  N 

N  N 

1 

II 

S  'LyR^yWyT’ 

y^/”/  _  ^ 

p^nj q^nj 

p^nj q^nj 

So  the  2-port  admittance  matrix  representing  the  jimction  is 


'i?' 

'^Joo  1 

YjOl 

-yf 

fL 

Yjll_ 

.v3. 
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where 

N  N 

^Joo  -  ~  L  L^L  Unv^L 

P*nj  q^nj 

N  N 

^J0\-~  L  ^inv^R 

p^nj  q^nj 

N  N 

y  _  ^  ^  Y^jPyP^Y^^^ 

^J\Q-~  Zj  Zj^R  Unv^L 

P*nj  q*nj 

N  N 

V  -  V  V  yf^sPYP^Y^^J 

^J\\-  ^R  “2-  Z^^R  Unv^R 

P*nj q*nj 

If  port  nj  of  Yl  and  Yr  are  connected  together  and  for  there  to  exist  a  solution  in  the  absence 
of  an  excitation  source  corresponding  to  a  mode  of  propagation,  then 

ll'  =-‘r  ‘"“‘''l  -K' 

Therefore 

Hence,  a  physical  value  for  is  obtained  when  the  following  transverse  resonance 
condition  is  satisfied 

t^/00  ^-/oi  +  ^^10  ^  =  0 

2. 3. 2.7  Numerical  Techniques  for  Solving  the  Trairsverse  Resonance  Condition 

To  complete  the  section  analysis,  the  propagation  constants  for  the  z-directed  modes,  or 
z-modes,  must  be  determined.  This  is  achieved  by  performing  an  iterative  search  for  kz 
values  which  satisfy  the  transverse  resonance  condition  derived  in  Section  2. 3. 2. 6.  Two 
markedly  different  search  methods  were  used  here  to  determine  the  k^  values  satisfying 
the  transverse  resonance  condition,  depending  upon  whether  the  required  z-modes  were 
propagating  or  evanescent  (purely  real  or  imaginary  kz),  or  complex  (complex  kz). 

Propagating  or  evanescent  modes 

To  find  the  propagation  constants  for  purely  propagating  and  evanescent  z-modes,  an 
iterative  search  for  both  zeroes  in  the  determinant  and  minima  in  the  minimum  singular 
value  of  [Y®"^]  was  conducted  by  decrementing  kz  from  some  initial  value.  The  minimum 
singular  value  of  was  found  using  matrix  singular  value  decomposition  [29],  a 

technique  first  applied  to  the  analysis  of  two-dimensional  microwave  and  millimetre- 
wave  structures  by  Labay  and  Bornemann  [30].  Figure  2.9  shows  a  typical  z-mode  solution 
corresponding  to  a  minimum  (zero)  in  the  minimum  singular  value  and  a  zero  in  the 
determinant  of  [Y*""'].  Note  the  determinant  pole  located  close  to  the  zero  crossing.  This 
was  a  common  feature  of  the  behaviour  of  the  determinant  near  z-mode  solutions,  and 
prompted  the  use  of  the  more  well-behaved  minimum  singular  value.  It  must  be  noted  that 
the  minimum  singular  value  was  not  entirely  reliable  on  its  own  because  of  extremely 
narrow  minima  sometimes  encountered.  Such  minima  cannot  be  detected  unless 
prohibitively  small  kz  search  step  sizes  are  used.  For  this  reason  both  the  determinant 
and  the  minimum  singular  value  of  [Y*^'”]  were  used  to  ensure  the  reliable  identification  of 
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valid  kz  solutions.  Instances  where  both  the  determinant  and  the  minimum  singular  value 
were  not  well-behaved  near  a  given  z-mode  solution  were  rare. 

The  following  strategy  was  adopted  for  finding  propagating  and  evanescent  z-modes: 


(i)  Start  with  a  value  of  no  larger  than  k^JEr  ,  where  is  the  maximum 

dielectric  constant  in  the  section.  Do  a  coarse  search  by  stepping  in  the  direction  of 
decreasing  kz  with  a  specified  step  size.  Look  for  minima  in  the  minimum  singular  value  or 
a  change  of  sign  in  the  determinant  of  Closely-spaced  zero-pole  pairs  in  the 

determinant  which  may  not  result  in  a  change  of  sign  over  a  coarse  step  can  also  be 
detected. 

(ii)  If  a  possible  solution  is  detected  during  the  coarse  search,  a  prediction  of  the 
location  of  the  solution  is  made  and  a  fine  search  initiated  over  a  range  of  one  coarse  step 
either  side  of  the  estimated  position  of  the  kz  solution.  The  fine  step  size  is  typically  set 
to  1/10  of  a  coarse  step.  As  for  the  coarse  search,  look  for  minima  in  the  minimum  singular 
value  or  a  change  of  sign  in  the  determinant  of 
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(iii)  If  a  possible  solution  is  detected  during  the  fine  search,  an  attempt  is  made  to  use  a 

quadratic  fit  to  the  minimum  singular  value  of  to  estimate  the  position  of  the 

solution.  If  this  estimate  corresponds  to  a  minimum  in  the  minimum  singular  value  then 
an  iterative  procedure  based  on  the  Brent's  method  is  used  to  refine  the  value  [31]. 
Otherwise,  a  quadratic  fit  is  performed  on  the  determinant  to  estimate  the  position  of  a 
possible  zero  crossing.  If  a  determinant  zero  crossing  is  detected,  then  the  estimate 
obtained  is  refined  using  a  procedure  based  on  the  method  of  false  position  [32]. 

(iv)  Before  acceptance  of  the  k^  value,  a  final  check  is  made  to  ensure  that  it 
corresponds  to  both  a  zero  crossing  in  the  determinant  and  a  minimum  in  the  minimum 
singular  value  of  [Y®""*].  If  the  solution  is  valid,  the  next  z-mode  is  sought  using  the  same 
procedure.  If  the  fine  search  is  not  yet  at  the  end  of  its  range  (one  coarse  step  beyond  the 
original  estimate  of  the  solution)  then  it  is  continued  to  this  limit  in  an  attempt  to  detect 
other  solutions  close  to  the  last  solution.  This  situation  occurs  often  in  structures  with 
almost  degenerate  z-modes. 

For  uniform  sections,  values  of  kg  corresponding  to  the  z-mode  cutoff  frequencies  can  be 
found  rather  than  the  z-mode  propagation  constants  at  a  given  frequency  as  described 
above  for  general  structures.  The  search  algorithm  was  identical  to  the  k2  search.  For  the 

m*^^  z-mode,  the  cutoff  frequency/^  was  used  to  calculate  the  propagation  constant  k  ^  at  a 
given  frequency  using  the  simple  relation 


where 


k()  — 


27tf 


and  kr, 


iTrfc 
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Since  the  z-mode  cutoff  frequencies  were  independent  of  frequency  for  uniform  sections, 
determining  the  cutoff  frequencies  or  kz  values  at  one  frequency  was  sufficient  to  analyse  a 
number  of  frequencies.  The  kz  values  obtained  at  the  first  frequency  were  readily  scaled  for 
any  number  of  different  frequencies,  so  that  a  z-mode  search  was  not  required  for  more  than 
one  frequency.  This  substantially  improved  computational  efficiency  when  analysing 
broadband  structures. 


Complex  modes 

For  finding  complex  z-modes,  ie.,  z-modes  with  complex  k  2  values,  a  contour  integration 
technique  based  on  Cauchy's  Theorem  is  applied  in  the  complex  kz  plane  to  search  for 
zeroes  of  the  determinant  of  [Y®“"^].  This  technique  consists  of  an  integration  around 
overlapping  circular  contours  in  the  complex  kz  plane  to  determine  the  number  of  zeroes  and 
their  location  within  the  contour.  This  procedure  was  initially  developed  for  solving 
matrix  equations  by  Delves  and  Lyness  [33],  and  applied  to  the  analysis  of  microwave 
structures  by  Sorrentino  and  Lampariello  [34]. 

Cauchy's  theorem  states  that 
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where  Nz  and  Np  are  the  number  of  zeroes  and  poles  of  the  function  f(z)  of  the  complex 
variable  z  within  the  contour.  Zj  is  the  position  of  the  i*^  zero  and  pi  is  the  position  of  the 
i*^  pole  within  the  contour.  Interpolation  using  a  polynomial  of  order  4  was  used  to 
estimate  the  derivative  f(z)  from  values  of  the  function  f(z).  By  setting  n=0,  Sq  =  Nz  -  Np 
is  determined.  Provided  there  are  not  too  many  zeroes  or  poles  inside  the  contour,  then 
with  a  sufficiently  high  order  of  n,  the  various  Zi  and  pi  can  be  determined  analytically. 
Appendix  A.5  describes  analytic  techniques  for  extracting  the  zero  and  pole  locations  from 
the  moments  Sn  for  a  limited  number  of  poles  and  zeroes  within  the  contour.  Larger 
numbers  of  poles  and/or  zeroes  were  not  considered  because  of  numerical  accuracy  problems 
which  stemmed  from  the  fact  that,  in  practice,  most  of  the  zeroes  and  poles  occurred  in 

closely-spaced  zero-pole  pairs,  so  that  z^  -  p^  was  small  for  n  =  1,  2,  3,  ....  This  was  a 

fundamental  limitation  on  the  effectiveness  of  this  technique.  However,  complex  modes 
were  generally  few  in  number  and  widely-spaced  in  the  kz  plane,  allowing  this  method  to 
be  used  effectively.  If  too  many  zeroes  or  poles  were  inside  a  given  contour,  smaller 
contours  could  be  used  to  cover  the  same  area  in  the  complex  plane. 

2. 3.2. 8  Section  Voltage  and  Current  Calculation 

Once  the  values  of  k^  are  found  for  each  z-mode  for  m  =  0  to  Nz,  the  voltages  and  currents 
on  the  equivalent  transmission  line  network  representing  propagation  in  the  x  direction 
can  be  determined.  These  can  then  be  used  to  calculate  the  actual  electric  and  magnetic 
fields  in  the  cross-section  for  each  z-mode,  and  hence  the  z-mode  coupling  at  section 
interfaces.  This  in  turn  allowed  sections  to  be  cascaded  together  to  construct  the  complete 
waveguide  element. 

The  first  step  is  to  use  the  transverse  resonance  condition  to  solve  for  the  modal  voltages  at 
the  junction 


This  equation  is  sufficient  to  solve  for  to  within  an  arbitrary  factor.  This  is  accounted 
for  by  allowing  one  of  the  modal  voltages  to  be  set  to  1.  If  voltage  Vp  is  set  to  1  and 
removed  from  the  unknown  voltage  vector,  then  the  p^^  coluiim  of  can  be  taken  to 

the  right-hand  side,  and  the  resulting  equation  solved  for  the  other  voltages  using  a  least 
squares  technique.  The  solution  is  ideally  unique,  but  the  least  squares  approach  ensures 
robustness  to  numerical  error.  In  matrix  terms,  the  equation  to  be  solved  is 


■^0(p-l)Yo(p+l)' 


.^sum  .^um 
N(p-l) ^  N(p+1)' 


If  this  equation  is  written  as  A.V  =  b,  then  the  least  squares  solution  would  be  found  by 
solving  A'*^A.V  =  A'*‘.b  using  standard  techniques  for  linear  matrix  equatiorrs  [35]. 


The  modal  voltage  matrices  at  the  other  ports  in  parallel  with  the  junction,  (n^f^nj),  are 
found  from  the  junction  voltages  using  this  equation  from  Section  2.3.2.6 
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P*nj 

where  N  is  the  number  of  ports  in  the  junction  plane  and 

The  currents  at  all  of  the  ports  in  parallel  with  the  junction  were  determined  using  Yl 

/£=  ii'2'’v'’=-4 

P*nj 

Once  the  voltages  and  currents  at  the  ports  in  parallel  with  the  junction  ports  were  found, 
the  voltages  and  currents  in  the  remainder  of  the  cross-section  could  be  determined.  For  a 
step-down  discontinuity  to  the  left  of  the  junction  plane,  the  voltage  coupling  could  be  used 
to  determine  the  voltage  on  the  left  side  of  the  discontinuity  from  the  known  voltage  on 
the  right  side.  The  current  to  the  left  of  the  discontinuity  was  calculated  by  evaluating 
the  total  admittance  of  that  portion  of  the  structure  from  the  left  side  boundary  up  to  (but 
not  including)  the  discontinuity,  Yieft/  so  that  I  =  Yieft  V  yielded  the  required  modal 
currents.  For  a  step-up  discontinuity  to  the  left  of  the  junction  plane,  the  current  coupling 
could  be  used  to  evaluate  the  modal  currents  to  the  left  of  the  discontinuity  from  the  known 
modal  currents  on  the  right  side.  As  for  the  step-up  case,  the  admittance  of  the  structure  to 
the  left  of  the  discontinuity,  Yieft/  was  calculated.  The  modal  voltages  on  the  left  side  of 

the  discontinuity  were  subsequently  determined  using  V  =  YleftI- 

Once  the  voltage  and  current  at  the  right  edge  of  a  region  Vright  ai'd  Iright  were  known,  it 
was  a  simple  matter  to  calculate  the  voltage  and  current  at  the  left  edge  Vieft  and  Iieft 
using 

^lefl  ~  Ylo[^right  “^ll^right] 

Iieft  =  YooVieft  +  YoiVright 

where  Y  is  the  subregion  admittance  matrix.  Note  that  Yjj  (i,j  =  0,1)  are  diagonal.  The 
next  discontinuity  (if  any)  is  located  at  the  left  edge  of  the  region,  so  this  procedure  is 
repeated  until  the  left  boundary  of  the  structure  is  reached  to  determine  the  j:-mode 
voltages  and  currents  for  the  portion  of  the  structure  to  the  left  of  the  junction  plane. 

The  modal  voltages  and  currents  to  the  right  of  the  junction  plane  are  calculated  in  a 
similar  fashion,  the  differences  being  that  the  voltage  coupling  was  used  to  calculate  the 
modal  voltages  across  step-down  discontinuities,  while  the  current  coupling  was  used  at 
step-up  transitions.  The  admittance  of  the  portion  of  the  structure  to  the  right  of  the 
discontinuity,  Yfight,  was  then  used  to  calculate  the  voltage  or  current  as  appropriate. 
Once  the  voltage  and  current  at  the  left  edge  of  a  subregion  was  known,  the  corresponding 
voltage  and  current  at  the  right  edge  were  determined  using 

Y right  =  Yq  1  [lleft  “  Yqq V ] 

Ifight  =  Yio  Yieft  +  YiiVright 

This  procedure  was  repeated  until  the  right  side  boundary  of  the  structure  was  reached. 
The  s:-mode  voltages  and  currents  across  the  entire  two-dimensional  cross-section  were  then 
known. 
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2. 3.2. 9  Section  Field  Calculation 

The  modal  voltages  and  currents  were  used  to  calculate  the  fields  in  the  cross-section,  as 
well  as  the  total  power  flowing  through  it.  The  y  and  z  components  of  the  fields  have  been 
derived  previously  in  order  to  evaluate  the  coupling  at  region  interfaces.  To  complete  the 
field  calculations,  the  x  component  of  the  electric  and  magnetic  fields  of  the  m*^  z-mode 

were  derived.  E™  was  determined  from  and  and  was  determined  from  and 

E^  using  Maxwell's  equations,  with  the  TE  and  TM  field  components  separated  as  follows 


jcoEQe^E^  =  V  X  H™  • 
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The  total  power  flowing  through  the  cross-section  in  the  m*^^  z-mode  was  determined  by 
integrating  the  z-directed  power  density  over  each  region  and  summing  over  all  the  regions 
in  the  cross-section 


P"  =  /  j  -  E^Hf)dy  & 

yi 

All  three  field  components  are  listed  in  Appendix  VI  for  TE  and  TM-to-x,  y,  and  z  modes, 
together  with  the  z-directed  power.  For  each  z-mode,  and  were  normalised  by 

setting  to  1  for  propagating  modes  (k2  real),  and  j  =  for  evanescent  modes  (kz 
imaginary)  and  complex  modes  (kz  complex). 

2.3.3  Three  Dimensional  (Element)  Analysis 

There  have  been  a  number  of  mode  matching  methods  developed  to  treat 
three-dimensional  discontinuities.  Early  work  concentrated  on  discontinuities  between 
rectangular  waveguide  sections  [36]-[42],  while  later  efforts  have  been  directed  toward 
more  complex  structures  based  on  firmed  waveguides,  cross-irises,  or  T-septa  [43]- [50]. 
Generally,  analysis  has  been  restricted  to  relatively  simple,  homogeneous  structures.  The 
more  general  derivation  presented  in  this  section  seeks  to  extend  the  capability  of  the 
mode  matching  technique  beyond  these  limitations  to  allow  accurate  modelling  of 
microstrip,  finline,  and  dielectric  waveguide  discontinuities  with  a  single,  generic 
analysis. 

The  two-dimensional  section  analysis  described  in  Section  2.3.2  yields  the  z-mode 
propagation  constants  and  modal  fields  for  a  given  number  of  z-modes  in  each  section  of  the 
structure  to  be  analysed.  It  remains  only  to  cascade  these  sections  together  (in  the  z- 
direction)  to  represent  the  entire  waveguide  element.  To  perform  this  cascading,  it  is 
necessary  to  first  calculate  the  coupling  between  z-modes  in  adjacent  sections  using  the 
modal  fields  calculated  according  to  the  expressions  in  Section  2.3.2.9.  Combining  this 
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coupling  with  the  scattering  matrices  for  each  section  allows  the  scattering  parameters  of 
the  entire  element  to  be  determined. 

2.3.3.1  Mode  Coupling 

The  z  modes  in  adjacent  sections  are  coupled  at  the  interface  between  the  sections.  This 
interface  represents  a  discontinuity  in  the  section  boundaries  or  dielectric  constants.  The 
z-mode  coupling  can  be  determined  from  continuity  of  the  tangential  electric  and  magnetic 
fields  in  the  {x-y)  plane  of  the  discontinuity  and  from  the  boundary  conditions.  Consider  a 
step  change  between  two  sections  A  and  B  of  an  element  as  shown  in  Figure  2.10. 


Using  the  two-dimerrsional  analysis  discussed  in  Section  2.2,  an  equivalent  z-mode  voltage 
and  current  can  be  defined  such  that  the  fields  tangential  to  the  z-direction  in  sections  A 
and  B  are  given  by 

m 

N* 

m 

n 

n 

where  the  summation  is  over  both  TE  and  TM  eigenmodes  and  propagation  is  assumed  to  be 
in  th  z-direction. 

Assuming  that  the  boundary  of  section  B  entirely  includes  the  boundary  of  section  A  at  the 
interface,  the  field  matching  at  the  discontinuity  requires  that 

^  ,  over  aperture 

~  LO,  otherwise 
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and 


Hr,  =H,^  over  aperture 

By  equating  the  modal  expansions  for  the  fields  in  regions  A  and  B  in  the  aperture,  the 

coupling  relations  between  V^^(z)  and  (z),  and  between  and  can  be 

determined.  This  procedure  was  followed  using  each  of  the  three  mode  types,  ie.,  TE, 
TM-to-JC,  y,  and  z.  The  coupling  expressions  were  formulated  in  such  a  marmer  as  to 
guarantee  conservation  of  complex  power  across  the  discontinuity  for  a  finite  number  of 
z-modes  in  each  section.  In  this  section,  the  z-mode  coupling  using  TE  and  TM-to-y 
eigenmodes  is  derived.  The  corresponding  derivation  for  TE  and  TM-to-x  and  TE  and 
TM  to-z  modes  is  given  in  Appendix  A.7. 
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Both  sides  of  Equation  2.1  were  multiplied  by 
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Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the 
y-modes  yields 


where 
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Both  sides  of  Equation  2.2  were  multiplied  by 
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Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the 
y-modes  yields 
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In  addition,  multiplying  both  sides  of  Equation  2.2  by 
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and  integrating  over  the  aperture  yields 
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Both  sides  of  Equation  2.3  were  multiplied  by 
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Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the 
y-modes  yields 
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y-component 


n!  Nf 


^  ^  .  mn  *n  ^  *  pq  *q 

I/I  I/a,  W''A,W  =  I/i  I/fi,  U)/«,W 


m  n 


Both  sides  of  Equation  2.4  were  multiplied  by 


j^TE  r  * 

k  0}jjQ  ,rs  *  ,s  * 

I^^/A,a)  v'a/}’) 


Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the 
y-modes  yields 


cv«,=Dv.^ 


where 


um  * 


A  K  (DfXQ  p  <np  <mp  n, 

Cmn  =  £  ^-72 - j  ^ 


P  xi 

.rTE  t^TE  „  * 

N  A  Nd  h-ffl  /■.■Ml 


A  Bjc  COPQ  f  <nq  <mp  f  'P  * 

Cm/.  =  1 1  I  /s,  U)/a,  (^)  *  j  v'b,  (y'lVA,  (y)  dy 


P  9  '"“Pa 

In  addition,  multiplying  both  sides  of  Equation  2.4  by 


__  „rs  *  i.J 

I/a,(^)  VA,(y) 


and  integrating  over  the  aperture  yields 


where 


MV«  =NVa 


fjTM  \jTE  Y 


=  ^  ^  (x)  dx\VA^{y)VAy{y)  dx 


P  9  Xi  yi 

^  ,nq  „mp  *  .<7  .p  * 

^mn  =  lh,(x)h,  (x)  dx\VB^(y)VB^iy)  dy 

P  9  Xi  yi 

Complex  Power  Conservation 

Several  relations  now  exist  between  the  z-mode  voltages  and  currents  at  section  interfaces. 
Ultimately,  a  single  coupling  expression  for  the  voltages  and  a  single  coupling  expression 
for  the  currents  is  required.  To  determine  the  means  by  which  the  previously  derived 
coupling  expressions  can  be  combined,  a  conservation  of  complex  power  technique  is  used. 
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Conservation  of  the  complex  power  flowing  in  the  z-direction  normal  to  the  section 
interface  is  enforced  regardless  of  the  number  of  z-modes  used. 

The  total  power  carried  by  the  structure  is  summed  over  all  z-modes  as  follows 

Xuyu 


^tJl 

^  m 


m  n 


m  n  xiyi 

Using  the  field  relations  in  Section  2.3.2.9,  the  power  flow  through  any  given  section  is 
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Therefore  for  section  A  and  section  B  can  be  expressed  in  terms  of  the  previously 
derived  coupling  matrices.  In  addition,  the  coupling  between  the  modal  fields  in  sections 
A  and  B  is  required  to  be  such  that  the  power  coupled  from  z-modes  in  section  A  to  z-modes 
in  section  B  is  the  same  as  the  power  coupled  from  z-modes  in  section  B  to  z-modes  in 
section  A. 


itot 


t  r  1 

= 

k: 

P  1 
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and 
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So  for  conservation  of  complex  power,  it  is  required  that 

D^+N^-L^=B  +  J-Q 
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This  relation  can  be  readily  verified  by  inspection  of  the  individual  terms  in  each  matrix 
listed  previously.  To  obtain  a  single  coupling  expression  for  the  z-mode  current,  and  a 
single  coupling  expression  for  the  z-mode  voltage,  write 

«Aj  =  [Cl][lBj.K]  =  [Cv]K^ 

Therefore,  for  power  conservation 

[Cif  [d  -h  M'*'  -  K'*']  =  [A  -1-  H  -  P]  [Cv] 

To  ensure  complex  power  conservation  regardless  of  the  number  of  z-modes  used  in  each 
section,  the  following  expressions  for  the  current  and  voltage  coupling  will  be  used 

[Ci]  =  [C-hM-Kr^  [D  +  N-L] 

[Cv]  =  [A  +  H-Pr^  [B  +  J-Q] 

=  [A-l-H-Pr^  [dUn^-L'^ 

The  corresponding  derivations  for  TE  and  TM-to-x  and  TE  and  TM-to-z  modes  are  presented 
in  Appendix  VII.  Appendix  VIII  contains  a  summary  of  the  coupling  expressions  for  each  of 
the  three  mode  options. 

For  sections,  with  uniform  dielectric  constants  the  self-coupling  matrices  [A  +  H  -  P]  and 
[C  +  M  -  K]  are  diagonal.  Therefore,  calculating  these  matrices  and  their  inverses  requires 
only  a  small  computational  effort  for  this  case. 

2.3.3.2  Section  Discontinuity  Scattering  Matrix 

Once  the  z-mode  voltage  and  current  coupling  matrices  are  determined  for  each 
two-dimensional  section  discontinuity,  the  scattering  matrix  of  the  discontinuities  can  be 
determined.  The  coupling  expressions  for  the  z-mode  voltages  and  currents  were  used  in  the 
form  used  previously  in  Section  2.3.2.4 


for  "step-up"  transitions  (section  A  cross  section  enclosed  by  section  B  cross  section),  and 


for  "step-down"  transitions  (section  B  cross  section  enclosed  by  section  A  cross  section).  Pj, 
and  Py  are  the  current  and  voltage  couplings,  respectively.  These  coupling  expressions  are 
listed  in  Appendix  VIII  for  the  various  choices  of  eigenmodes. 


As  for  the  two-dimensional  case  discussed  in  Section  2.3.2.4,  note  that  the  currents  in  the 
smaller  section  could  be  expressed  in  terms  of  the  currents  in  the  larger  section,  but  not  vice 
versa.  This  is  because  there  was  no  boundary  condition  for  the  magnetic  field  on  the  wall  of 
the  aperture.  Similarly,  the  voltage  in  the  larger  region  could  be  determined  from  the 
voltage  in  the  smaller  region  and  the  boundary  constraint  that  the  component  of  the 
electric  field  tangential  to  the  section  interface  must  disappear  on  the  aperture  wall. 
However,  as  for  the  two-dimensional  case,  this  boundary  constraint  means  that  a  coupling 
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expression  relating  the  voltage  in  the  smaller  region  to  an  arbitrary  voltage  in  the  larger 
region  carmot  be  defined,  since  an  arbitrary  electric  field  cannot  be  specified  in  the  larger 
region. 

A  scattering  matrix  relating  the  amplitudes  of  the  incident  and  reflected  waves  on  each 
side  of  the  discontinuity  is  required  and  this  can  be  derived  from  the  voltage  and  current 
coupling  relationships.  This  scattering  matrix  is  written  as 


section  interface 


Figure  2.11  Single  input  section  coupled  to  N  output  sections 
Step-down  transitions 

Figure  2.11  schematically  shows  a  single  input  section  coupled  to  N  output  sections  through 
a  step-down  discontinuity.  The  voltage  and  current  coupling  between  the  input  section  and 

the  n*^  output  section  are  written  as  Py  ^  and  respectively.  Representing  the  incident 
and  reflected  wave  amplitudes  as  shown  in  Figure  2.11,  the  following  relations  can  be 
defined 

n 
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Therefore 


b'  =  +  +  b" -c‘ 


This  leads  to 


(>'  =  -  U  +  X  Py ’P]'’*  U  -  c' 


+2  U  +  SPWP</* 


In  matrix  form 


^//^pW  U+  U  +  ^pWp(p)  U-I^Pif’pW  c' 


+1;  ua^-2P<">  u+|;p</>pi'’>  Pi?*  b'J 


[Sill  I  [  S12  ]  C 


S21 


where 


Sii  =  - U  +  XPv'”P;'’’  U-2pV'’’P;'’* 


Si2„=2  U  +  Xp'/’P/*”  Pi?* 


S2i„  =  Pi"*[U-Siil 
S22„„=[U5„„-Pi“*Sl2„l 
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Step-up  transitions 

Figure  2.12  schematically  shows  N  input  sections  coupled  to  a  single  output  section  through 
a  step-up  discontinuity.  The  voltage  and  current  coupling  between  the  n*^^  input  section  and 

the  output  section  are  written  as  Py  ^  and  respectively.  Representing  the  incident  and 
reflected  wave  amplitudes  as  shown  in  Figure  2.12,  the  following  relations  can  be  defined 

Therefore 


section  interface 


Figure  2.12  N  input  sections  coupled  to  a  single  output  section 
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In  matrix  form 


aJ  =  i  -2P<">  U  +  XPif>Pi'’*  pS,’’  4 


+p5''*  U+  U  +  £p</’p{'’'  U-^P^/'p*/’  b" 


where 


[  S21  ] 


[U5mn-Pf^S2lJ 

pf'^[U-S22] 


=  2  U  + 


’)  oW 


N  1-1  r  N 

-  u  +  X  u-X  P^^^Pf"^ 

p  J  L  p 


2.3.3.3  Scattering  Matrix  Cascade 

The  modal  scattering  matrices  for  each  two-dimensional  section  discontinuity  (in  the  x-y 
plane)  can  be  cascaded  in  the  z-direction  with  the  simple  2-port  scattering  matrices  of 
each  section  to  produce  the  overall  scattering  matrix  for  the  entire  waveguide  element. 
The  general  case  of  two  cascaded  multiport  scattering  matrices  is  shown  in  Figure  2.13. 

The  incident  and  reflected  wave  amplitudes  (cl,bl),  (cll,bll),  and  (clll,blll)  are  related  by 

(b']j^r(s[i]  [Si2iT[p'J 
[c"]  [Sill  [s^2m*>"i 


[b"l  _r[sf'|]  [sfj]  [c"] 


[<^"'1  N']  [S?2l  [b'"] 
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Figure  2.13  Multiport  scattering  matrix  cascade 

The  required  cascaded  scattering  matrix  is  given  by 

[b'll  r[sfi)  [sf2]T[c'] 

[[S5i]  [S|2]|[b“l 

The  parameters  at  the  cascaded  ports  are  eliminated  by  first  expressing  [b^^]  as  follows 
[b'']  =  [Sf,][c"]  +  [S"][b'"] 

=[s;!l[s;,][c']qs^][b'']]+[s.»][b»'] 

and  then  by  eliminating  [c^^]  to  obtain 

[b"]=[[s;',]+[s;,](u-[s"][s;,])‘'[s;;][s;,]][c'] 

+K](u-[s;;][s',])'‘[s;'|b"'] 

and  finally  by  writing  an  expression  for  as 

[c"'l  =  [S"|][c"]  +  [S"][b'l 

=[s",][[s;,][c"]+[s;,][b"j]+[s"][b'] 

= [s",  ][u + [s;  j(u  -  [s;;  is',  if'is"  i][b"'is;,  ic'i 

+[[s"]+[s"isy(u-[s;;is;,])"‘[s"]][b"'] 
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the  required  parameters  can  be  identified  as 


Similar  expressions  for  the  case  of  two-port  scattering  matrix  cascading  are  presented 
in  [37],  The  treatment  shown  above  is  superior  to  that  described  in  [37]  in  that  multiport 
scattering  matrices  can  be  handled,  and  onJy  a  single  inverse  needs  to  be  evaluated. 


The  2-mode  scattering  matrices  for  uniform  sections  are  diagonal,  with  elements 


Qtnn  _  amn  _  n 
>^11  -*^21 

Sn  =  ^21  =  Dmn  =  Smn  exp(-7fcf  4) 

where  Ig  is  the  length  of  the  section  and  bmn  is  the  Kronecker  delta.  Note  that  this  is  the 
only  point  where  the  section  lengths  appear,  since  all  other  quantities  relating  to 
propagation  constants  and  mode  coupling  are  related  to  the  two-dimensional  cross  sections 
of  the  various  sections  in  the  element.  This  means  that  only  the  final  scattering  matrix 
cascade  needs  to  be  recalculated  if  section  lengths  are  varied.  This  leads  to  considerable 
savings  in  computational  effort  where  optimisation  of  a  structure  is  required  with  respect 
to  the  lengths  of  individual  sections. 

To  perform  the  section  cascade,  the  two-port  scattering  matrices  for  the  uniform  sections 
are  first  cascaded  with  the  scattering  matrix  of  the  discontinuity  immediately 
preceding  each  section.  This  yields  a  multiport  scattering  matrix  which  includes  the 
discontinuity  and  the  uniform  section  following  it  which  is  given  by 

[s;'i]=[sfi] 

(S('2]  =  [Sf2)[D] 

[Sf,]  =  [D][S?,] 

[Sf2l  =  [D][S?2][Dl 

As  described  in  [37],  strongly  evanescent  modes  for  which  I  /g  I  <  15,  where  Ig  is  the 
length  of  the  section,  are  not  considered  in  the  cascade.  While  such  modes  play  an 
important  role  in  evaluating  the  z-mode  scattering  at  individual  discontinuities,  their 
fields  are  extremely  weak  away  from  discontinuities.  Therefore,  once  they  have  been 
accounted  for  in  the  coupling  evaluation,  their  subsequent  effect  on  the  overall  scattering 
matrix  cascade  is  negligible. 

The  section  cascade  bemns  at  the  specified  input  end  of  the  element,  with  successive 
sections  represented  by  ^  cascaded  to  the  end  of  the  previous  section.  The  total  cascade  up 
to  the  section  to  be  added  is  represented  by  S^.  The  general  multiport  cascade  is  thus 
applied  repetitively  to  construct  the  entire  element.  The  input  and  output  ports  for  the 
element  are  assumed  to  be  terminated  with  matched  impedances.  These  scattering 
parameters  constitute  the  output  of  the  three-dimensional  analysis. 
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3  RESULTS 


3.1  Two-dimensional  Structures 

Various  two-dimensional  waveguiding  structures  for  which  theoretical  or  experimental  data  are 
known  were  analysed  using  the  transverse  resonance  technique  described  in  Section  2.  By 
comparing  the  results  of  the  analysis  with  existing  results,  both  the  transverse  resonance  theory 
and  its  implementation  in  the  program  TERESA  (for  TransversE  RESonance  Analysis)  could  be 
tested.  In  each  structure,  TE  and  TM-to-y  eigenmodes  were  generally  used,  although  TE  and 
TM-to-z  modes  are  more  appropriate  for  structures  with  a  uniform  dielectric  constant.  Although 
affecting  the  required  amount  of  computational  effort,  the  particular  choice  of  TE  and 
TM-to-x,  y,  or  z  eigenmodes  was  found  to  have  no  effect  on  the  results. 

3.1.1  Double  Ridged  Waveguide 

Montgomery  [9]  calculated  the  cutoff  frequencies  of  TE  modes  in  the  double  ridged 
waveguide  shown  in  Figure  3.1.  A  comparison  was  made  these  results  and  results  obtained 
from  the  transverse  resonance  analysis.  This  data  is  shown  in  Table  3.1.  In  the  side 
regions,  9  TE  and  8  TM-to-z  modes  were  used,  while  3  TE  and  2  TM-to-z  modes  were  used  in 
the  gap.  The  number  of  modes  in  the  side  regions  was  determined  from  relative 
convergence  given  the  number  of  modes  in  the  gap.  The  number  of  gap  modes  was  chosen 
according  to  convergence  of  the  final  results.  Note  that  one  more  TE-to-z  mode  was  used 
than  TM-to-z  modes.  This  is  because  the  lowest  order  TE-to-z  mode  is  the  TEM-to-u  mode 
for  which  ky  =  0.  This  mode  has  no  TM-to-z  counterpart,  so  that  one  extra  TE-to-z  mode 
must  be  used  to  attain  the  same  ky  value  as  the  highest  order  TM-to-z  mode.  This 
situation  also  arises  when  TE  and  TM-to-x  modes  are  used.  If  TE  and  TM-to-y  modes  are 
used,  one  more  TM  mode  than  the  number  of  TE  modes  is  used  to  account  for  the 
TEM-to-u  mode.  This  difference  arises  the  different  definitions  of  the  various  mode  types. 


Figure  3.1  Double  ridged  waveguide  cross  section  [9]  (dimensions  in  mm) 

The  transverse  wavenumbers  kx  for  various  modes  are  shown  in  Table  3.1.  The  modes 
labelled  as  "trough"  modes  are  so  named  because  the  fields  in  the  gap  are  only  weakly 
coupled  to  the  fields  in  the  side  regions.  The  agreement  in  Table  3.1  is  generally  better 
than  0.2  %.  The  simplicity  of  the  ridged  waveguide  cross  section  belies  its  importance  in 
many  applications  including  broadband  transitions  and  filters. 
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Table  3.1  Comparison  of  ridged  waveguide  TE  mode  cutoff  frequencies 


Mode 

f^^^^  (GHz) 

f™  (GHz) 

TEio 

6.857 

6.871 

TEio  (trough) 

15.104 

15.104 

TE20 

24.858 

24.890 

TE20  (trough) 

29.536 

29.541 

TE30 

32.024 

32.026 

TEii  (trough) 

33.271 

33.271 

3.1.2  Unilateral  Finline 

Finlines  are  an  important  structure  extensively  used  in  broadband  millimetre-wave 
applications  due  to  their  compactness  and  the  ease  of  mounting  active  devices.  The 
theoretical  results  of  Vahldieck  [21]  for  the  unilateral  finline  shown  in  Figure  3.2  were 
used  to  verify  the  transverse  resonance  results  for  this  guide.  This  structure  was  more 
complex  than  the  ridged  waveguide  in  that  it  contained  both  "step-up"  and  "step-down" 
discontinuities  on  either  side  of  the  slot  region,  as  well  as  a  dielectric.  Symmetry  about 
the  axis  shown  in  Figure  3.2  was  invoked  to  reduce  the  computational  effort  required.  For 
convergence,  2  TE  and  3  TM-to-y  modes  were  used  in  the  slot  region,  with  13  TE  and 
14  TM-to-y  modes  in  the  dielectric  region  and  8  TE  and  9  TM-to-y  modes  in  the  side  regions. 


Figure  3.2  Unilateral  finline  cross  section  [21]  (all  dimensions  in  mm) 

A  comparison  of  the  normalised  wavenumbers  kz/ko  is  shown  in  Figure  3.3.  Note  that  the 
agreement  is  generally  better  than  3  %.  The  discrepancies  between  the  two  sets  of  results, 
particularly  for  higher  order  z-modes,  may  be  due  to  numerical  inaccuracy  or  convergence 
difficulties  in  Vahldieck's  method  as  discussed  by  Mansour  and  MacPhie  [27].  Graph 
reading  errors  in  obtaining  Vahldieck's  results  could  also  contribute  to  the  observed 
discrepancy. 
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3.1.3  Suspended  stripline 

The  suspended  stripline  from  [27]  and  [51]  shown  in  Figure  3.4  was  analysed  using  the 
transverse  resonance  technique.  Symmetry  about  the  axis  shown  was  used  to  reduce  the 
required  computational  effort.  For  convergence,  3  TE  and  4  TM-to-y  modes  were  used  in  the 
narrow  region  imder  the  strip,  with  4  TE  and  5  TM-to-y  modes  in  the  other  regions. 


Figure  3.4  Suspended  stripUne  cross  section  [27] 
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The  results  were  compared  with  theoretical  results  published  by  Mansour  and  MacPhie 
[27],  and  Bornemann  [51]  as  shown  in  Figure  3.5.  The  agreement  with  [27]  is  excellent, 
whereas  some  discrepancies  are  again  noticed  with  Bomemarm's  results,  particularly  for 
the  higher-order  modes.  As  with  the  unilateral  finline,  numerical  inaccuracy  and 
convergence  difficulties  inherent  in  Bornemann's  method  are  the  most  probable  causes  of 
error. 

3.1.4  Coplanar  Waveguide 

The  coplanar  waveguide  is  another  structure  commonly  used  for  millimetre-wave 
applications.  The  coplanar  waveguide  shown  in  Figure  3.6  was  analysed  and  the  results 
compared  with  those  of  Mansour  and  MacPhie  [27],  and  Bomemann  and  Arndt  [13].  The 
calculation  of  the  dominant  and  first-order  odd  mode  normalised  propagation  constants 
kz/kfl  was  performed  with  the  depth,  d,  of  the  groove  set  to  0  and  0.5  mm.  Symmetry  was 
used  to  reduce  the  size  of  the  calculation,  with  an  electric  wall  inserted  on  the  axis  shown 
to  generate  the  desired  odd  symmetry.  For  convergence,  2  TE  and  3  TM-to-y  modes  were 
used  in  the  slots,  with  8  TE  and  9  TM-to-y  modes  in  the  side  regions.  For  d  =  0,  8  TE  and 
9  TM-to-y  modes  were  used  in  the  dielectric  regions,  and  13  TE  and  14  TM-to-y  modes  for 
the  d  =  0.5  case,  as  required  by  the  relative  convergence  criterion. 


For  both  groove  depths,  the  agreement  with  [27]  and  [13]  is  close,  as  shown  in  Figure  3.7, 
although  the  agreement  is  closer  for  the  d  =  0  case.  Note  that  increasing  the  groove  depth 
decreased  the  first-order  odd  mode  cutoff  frequency,  reducing  the  single-mode  bandwidth 
of  the  guide  [27]. 

3.1.5  Groove  Nonradiative  Dielectric  Guide 

Nonradiative  dielectric  waveguide  is  attractive  for  millimetre-wave  applications 
because  of  its  simplicity,  ease  of  fabrication,  and  low  loss.  Tong  and  Blundell  [52]  included 
grooves  in  the  upper  and  lower  plates  to  fix  the  dielectric  slab  in  place,  as  shown  in 
Figure  3.8.  The  inclusion  of  a  groove,  however,  renders  simple  single-mode  transverse 
resonance  modelling  inaccurate  because  of  mode  scattering  at  the  groove  discontinuities,  so 
that  a  generalised,  multimode  treatment  is  required. 
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Transverse  resonance  analysis  of  the  groove  nonradiative  dielectric  (GNRD)  guide  shown 
in  Figure  3.8  yielded  values  for  the  normalised  wavenumber  kz/ko  of  the  LSMn  mode 
which  could  be  compared  with  experimental  data  published  by  Tong  and  Blundell.  GNRD 
guides  with  relative  dielectric  constants,  Er,  of  2.56  and  4.0  were  analysed.  TE  and  TM-to-x 
modes  were  used  in  the  calculation,  with  5  TE  and  4  TM  modes  used  in  each  region  to  ensure 
convergence.  The  agreement  with  the  experimental  data  shown  in  Figure  3.11  is  generally 
better  than  5  %.  This  disagreement  is  also  evident  in  Tong  and  Blundell’s  calculations,  and 
so  probably  originates  from  experimental  errors. 


Figure  3.9 


Comparison  of 


kz 

ko 


versus  frequency  results  for  GNRD  waveguide 


.6  Shielded  Dielectric  Image  Guide 

A  comparison  was  made  between  the  theoretical  results  presented  by  Strube  and  Arndt  [40] 
and  the  results  from  the  transverse  resonance  analysis  for  the  shielded  dielectric  image 
guide  shown  in  Figure  3.10.  This  guide  differs  from  those  discussed  previously  in  that  it 
possesses  a  non-uniform  region.  For  this  reason,  TE  and  TM-to-y  modes  must  be  used.  For 
convergence,  5  TE  and  6  TM-modes  were  required. 


Figure  3.10  Shielded  dielectric  image  guide  cross  section 
(all  dimensions  in  mm) 
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Unlike  the  structures  previously  examined,  this  guide  was  found  to  support  a  complex  mode 
between  13.8  and  16.3  GHz,  as  can  be  seen  from  Figure  3.11.  Complex  modes  occur  as  a  result 
of  coupling  between  degenerate  z-modes.  Although  only  one  mode  is  shown  in  Figure  3.11, 
they  must  exist  in  degenerate  pairs  with  propagation  constants  of  opposite  sign 
(±  kz  =  ±  P  ±  j  a).  While  a  real  power  flow  exists,  the  total  power  transmitted  through 
the  cross-section  is  zero,  since  equal  amounts  of  power  are  transmitted  in  the  forward  and 
backward  directions.  Hence,  complex  modes  do  not  correspond  to  power  dissipation,  and 
are  supported  in  the  lossless  waveguide  model.  The  agreement  between  the  two  sets  of 
results  shown  in  Figure  3.11  is  excellent. 


3.2  Three-dimensional  Structures 

The  calculation  of  the  propagation  constants,  kz,  and  the  modal  fields  for  the  z-modes  in 
two-dimensional  sections  allowed  the  scattering  parameters  for  three-dimensional  waveguide 
elements  to  be  calculated.  As  for  Section  3.1,  structures  for  which  independent  theoretical  or 
experimental  results  could  be  found  were  analysed. 

3.2.1  Rectangular  Waveguide  Transformer 

The  simplest  type  of  element  to  analyse  comprises  cascaded  rectangular  waveguide 
elements.  A  typical  example  of  such  a  structure  is  the  Ku  to  X-band  transformer  shown  in 
Figure  3.12,  with  dimensions  as  shown.  This  structure  was  analysed  by  Arndt  et  al.  [39]. 
For  analysis  by  the  TERESA  program,  6  TE  and  5  TM-to-z  y-modes  with  60  z-modes  were 
used  in  each  section.  A  comparison  of  the  reflection  coefficient  I  Sn  I  given  in  [39]  and  by 
the  TERESA  program  is  shown  in  Figure  3.13.  The  small  error  between  the  two  sets  of 
results  can  be  attributed  to  reading  error  in  obtaining  the  results  from  [39]  graphically. 
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Figure  3.12  Four-section  Ku  to  X-band  waveguide  transformer  (all  dimensions  in  nun) 


3.2.2  Rectangular  Waveguide  E-Plane  Stubs 

A  comparison  of  the  magnitudes  of  Sn  and  S12  and  the  phase  of  Sn  with  known  results 
was  made  for  the  rectangular  waveguide  with  E-plane  steps  shown  in  Figure  3.14.  This 
structure  was  independently  analysed  by  Rozzi  and  Mongiardo  [53].  The  results  from  [53] 
are  compared  with  the  TERESA  output  in  Figure  3.15.  The  two  sets  of  results  agree  closely, 
except  at  14.75  GHz  where  TERESA  appears  to  predict  a  resonance.  Note  however  that 
the  Rozzi  and  Mongiardo  technique  also  shows  some  anomalous  behaviour  near  this 
frequency.  For  the  TERESA  analysis,  4  TE  and  3  TM-to-z  modes  and  35  z-modes  were  used 
in  the  smaller  sections,  with  6  TE  and  5  TM-to-z  modes  and  55  z-modes  used  in  the  larger 
sections.  Symmetry  was  used  to  reduce  the  size  of  the  calculation. 
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3.2.3  E-Plane  Insert  Filter 

The  E-plane  metal  insert  filter  shown  in  Figure  3.16  was  analysed,  and  the  results 
compared  to  calculations  by  Vahldieck  et  al.  [44],  and  measurements  by  Tajima  and 
Sawayeuna  [54],  This  comparison  is  shown  in  Figure  3.17.  The  results  from  the  TERESA 
analysis  agree  well  with  those  from  [44].  A  small  discrepancy  is  apparent  with  the 
measured  results  of  [54]  which  can  be  attributed  to  graph  reading  error.  For  the  TERESA 
analysis,  6  TE  and  5  TM-to-z  modes  and  55  z-modes  were  used  in  each  section.  Symmetry 
about  the  y-z  plane  was  invoked  to  reduce  computational  effort. 
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Figure  3.16  Rectangular  waveguide  E-plane  metal  insert  filter  (all  dimensions  in  mm) 
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3.2.4  Square  Waveguide  Iris  Polariser 

The  square  waveguide  iris  polariser  shown  in  Figure  3.18  was  analysed  and  a  comparison 
made  with  independent  calculated  results  from  [55].  These  components  are  used  to  excite 
circularly  polarised  waves  in  square  aperture  antennas  by  producing  a  90  degree 
differential  phase  shift  between  the  orthogonal  TEio  and  TEqi  modes. 


Figure  3.18  Square  waveguide  iris  polariser 
(all  dimensions  in  mm) 


For  the  TERESA  analysis,  20  TE  and  21  TM-to-y  modes  and  60  z-modes  were  used  in  each 
section.  Horizontal  and  vertical  symmetry  were  used  to  reduce  the  size  of  the  calculation. 
For  the  TEio  mode,  a  horizontal  magnetic  wall  and  a  vertical  electric  wall  were  used,  and 
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a  horizontal  electric  wall  and  a  vertical  magnetic  wall  for  the  TEqi  mode.  Figure  3.19 
shows  the  good  agreement  between  the  TERESA  results  and  the  analysis  presented  in  [55]. 
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3.2.5  Ridge  Waveguide  Notch 

The  notched  ridge  waveguide  shown  in  Figure  3.20  was  analysed  with  the  results 
compared  to  experimental  data  obtained  by  the  present  authors.  A  resonance  was  observed 
at  6.55  GHz.  The  close  agreement  between  theory  and  experiment  is  shown  in  Figure  3.21. 
For  the  TERESA  analysis,  60  z-modes  were  used  in  each  section  to  ensure  convergence. 
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3.2.6  Cross-Iris  Filter 

The  transmission  coefficient  I S21 1  was  calculated  the  cross-iris  filter  shown  in  Figure  3.22 
using  the  TERESA  analysis.  In  the  rectangular  section,  8  TE  and  9  TM-to-y  modes  with 
60  z-modes  were  used,  For  the  irises,  3  TE  and  4  TM-to-y  modes  were  used  in  the  smaller 
side  regions,  and  5  TE  and  6  TM-to-y  modes  in  the  larger  centre  region,  with  a  total  of 
30  z-modes  required  for  convergence.  Symmetry  about  the  vertical  and  horizontal  planes 
was  used  to  reduce  the  size  of  the  calculation.  The  results  obtained  compared  well  with 
the  calculated  and  measured  results  determined  by  Ihmels  and  Arndt  [50],  as  shown  in 
Figure  3.23. 


4  CONCLUSIONS 

A  three-dimensional  mode  matching  analysis  which  uses  the  generalised  transverse  resonance 
technique  to  analyse  the  constituent  two-dimensional  cross-sections  in  passive  waveguide  structures 
has  been  developed  and  successfully  tested.  Comparison  with  a  broad  range  of  waveguide  components 
for  which  independent  experimental  or  theoretical  data  could  be  found  revealed  close  agreement  in  all 
cases. 

This  method  potentially  offers  advantages  in  speed  and  efficiency  over  more  generally  applicable,  but 
computationally  intensive,  methods  such  as  the  finite  element  or  finite  difference  methods.  It 
promises  to  be  especially  useful  for  the  analysis  and  optimisation  of  broadband  waveguide  components. 

Further  work  will  concentrate  on  the  analysis  of  more  complex  three-dimensional  waveguide 
components  for  the  design  of  filters  and  broadband  transitions. 
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APPENDIX  I  MODE  FUNCTION  GENERATOR 

The  field  in  each  subregion  of  a  given  two-dimensional  section  within  the  element  to  be  analysed  is 
represented  by  an  expansion  in  terms  of  y-modes  which  are  TE  and  TM  to  the  x,  y,  or  z  reference 
directions.  The  choice  of  which  direction  was  most  appropriate  for  a  particular  type  of  structure  is 
discussed  in  Section  2.3.2.  The  general  technique  used  to  generate  the  fields  corresponding  to  these 
mode  functions  follows  Felsen  and  Marcuvitz  [10],  and  is  derived  in  this  appendix  from  Maxwell's 
equations.  These  equations  are 


VxH 
VxE 
V  •  eE 
V«pH 

From  Equation  I.l 

VxVxH 
(V  X  V  X  H)  X  Uz 

Using  the  identity 

V(A»B) 

then  (with  A  =  E  and  B  =  u^) 

[V  X  (jcoeE)]  X  Uz 


jcoeE 

-joopH 

0 

0 


(1.1) 

(1.2) 
(1-3) 
(1.4) 


V  X  (jcoeE) 

[V  X  (jcoeE)]  X  Uz 


(A»V)  B  +  (B»V)  A  +  A  X  (VxB)  +  Bx  (VxA) 


jcoe  [Ex  (V  X  Uz)  +  (E«V)  Uz  +  (Uz*V)  E  -  V(E»Uz)] 
jcoe[0  +  0  +  ^+^-VtEz-VzEz] 

j(oe[^-VtEz]  (1.5) 


where  Vt  is  the  transverse  gradient  operator,  and  Et  is  the  field  in  the  plane  transverse  to  the  reference 
direction  to  which  the  modes  are  to  be  TE  or  TM.  Note  that  since  the  modes  are  to  be  recalculated  for 
each  subregion  layer  of  differing  dielectric  constant,  e  is  regarded  as  uniform.  From  Equations  I.l  and 
1.2 


where 

From  Equations  1.5  and  1.6 


Now  from  Equation  I.l 

Thus 


V  X  V  X  H  =  k^H 
k2=  co^pe 


3E, 

dz 


X  Uz  -F  VjE, 

jcoe 


Vtx  Ht  =  jcoeEz 


Vt  X  Ht  •  Uz  =  Vt  •  Ht  X  Uz  =  jcoeEz 
VtVt  •  Ht  X  Uz  =  jcoeVtEz 


(1.6) 


(1.7) 
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From  Equations  1.3  and  1.4 


1“  =  j|j[k=H,xu,4.V,V,.H,xuJ 


= 


lt+1^ 


•  Ht  X  Uz 


(1.9) 


where  It  is  the  transverse  unit  dyadic  vix^x  +  ^y^y- 
Similarly, 

aHt 

dz 


=  -i(0£ 


•  u^xEt 


(I.IO) 


Let 


Mx,y,z)  =  ^Vi(z)ei(x,y) 
i 

Ht(x,y,z)  -  ^Ii(z)hi{x,y) 


Substituting  in  Equations  1.9  and  1.10 


y  ^e,  (x,y)  =  -jCOfl 

.  dz 


'^^hi(x,y)  =  -jwel 
:  dz  ' 


1 

Ir  +  ^2 

VrV,' 


These  equations  will  be  satisfied  if 


1  dVj 
Ii(z)  dz 

1  d/i 


^i(x,y)  =  -jm 


h/(j:,y)  =  -j(oe\ 


.  i 

1  +2*- 


Vi(z)  dz 

Now  RHS  is  independent  of  z  coordinate  in  both  cases  and  thus  we  can  write 

1  dV; 


•  ^/,(z)ht(x,y)xu2 

i 

•  U2X^V,(z)e,(x,y) 

i 

•  hiix,y)xu^ 

•  Uj  xe/(j:,y) 


liiz)  dz 
1  dlj 
Viiz)  dz 


le. 


ei(x,y)  =  -jKiZiei(x,y) 

hiix,y)  =  -jKiYihi(x,y) 


^  =  -yV,Z,7;(z) 

dz 
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APPENDIX  II  MODE  FUNCTION  DERIVATIONS 

Using  the  mode  function  generator  taken  from  Felsen  and  Marcuvitz  [10]  and  derived  in  Appendix  I,  the 
actual  y-mode  fields  can  be  calculated.  These  modes  are  defined  as  TE  or  TM  to  a  reference  direction  x, 
y,  or  z,  depending  on  the  nature  of  the  section  being  analysed.  TE  and  TM-to-y  mode  fields  are  derived 
in  Section  2.3.2.1,  and  this  appendix  contains  exactly  the  same  derivation  for  TE  and  TM-to-x  and  TE 
and  TM-to-z  modes.  For  x  and  z  reference  directions,  note  that  the  single  prime  denotes  TM  modes,  and 
the  double  prime  denoted  TE  modes.  This  is  opposite  to  the  convention  adopted  for  the  y  reference 
direction  in  Section  2.3.2.I.  The  reason  for  this  was  to  maintain  as  many  similarities  as  possible 
between  the  derivations  in  this  appendix  and  the  derivation  in  Section  2.3.2. 1.  The  different 
convention  is  also  convenient  in  that  it  reflects  the  manner  in  which  the  TERESA  source  code  was 
written  to  implement  the  layer  mode  function  calculation  for  each  subregion  in  a  given  two-dimensional 
section.  Note  that  the  y-mode  normalisations  given  in  Section  2.3.1  are  assumed. 


TE-to-x  {tx=  0) 


kxYj^\i'l{y,z)  =  (oe\ 

Thus  (after  cancelling  ei^z^  factor) 

1 


eT(y,z)  =  ey(y)e  ^'^z^Uy  +  e^(y)e 
hxCj.z)  =  h^(y)e~j’^z^Uy  +  h^(y)e~'^’^z^Uj^ 

VtVt 


If  +  ■ 


T''T 


UxXexCy.z) 


kJxihyUy+h^u^)  =  - 


fr  ^2 


O)ll0 


By 


j  +  koEr  fz^y  +  jkz  "  (^0  )ey^z  “ 


By 


So 


kx^x^ 


(OHq 


[k^Er-ky^z'^  jk; 


"  ^ 


Be 


Now,  since  V.e"  =  0,  and  ^  =  0 


Therefore 


If  is  chosen  to  be 


then 


=  ;Vz 


O/XO 


tf 


M 


kx 

COflQ 
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Similarly 


"  "  if/?  9 \  "  9^ 


(OjlQ' 


) 


From  V.e  "=  0 


-k^e  -  ik 

y  y  ■  ^  ^  dy 


So 


kx^x^z  ~ 


COliQ 


Hence 


-  By 


As  for  TM-to-y  modes,  set 
Using 


By  =  I  iy)B  ^  z 


-jk^z 


dl 


^  =  -jkyY  V  (y)  =  -jweoV  (y) 
dy 


leads  to 


It  11 


J 


TM-to->:  (h;c=  0) 


e7’(>>,z)  =  e^(y)e  ^’^z^Uy  +  B^{y)B 


kj^Z^ej'(y,z)  =  CO 
Thus,(after  cancelling  e'i*^2^  factor) 


hT(y,z)  =  hy(y)B  •^Vuy +  /z^(y)e 

•  hjiy,z)'X 


cy  £(y) 


(0£Q£f 


{ko£Q -ky)h^+jk^^ 


Now,  since V.h'  =  0  and  hr  =  0 


By 


AK 


So 


^x^x^y  ~ 


Q)£Q£f 
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If  Xj.  is  chosen  to  be 


then 

Similarly 

From  V.h'  =  0 

So 

Hence 

As  for  TE-to-y  modes^  set 
Using 

leads  to 


03£q£j- 


1 


0}£Q£r  I 


]£Q-k'^'jhy+jk. 


dh^ 


.k^h  =  ik  ^ 

Kyriy  JK^ 


(oeoEr 

t  t 

G-y  —  "h,/ 


X  * 


9V 


V  (y)  e-jl^zz 

-jkyYI(y)  =  -jcoi^oUy) 

t  i  9h,, 


kz  3y 

“•^0  tV  ^  -ik  z 
-j—  I  (y)  e 


e'jkz^ 


TE-to-z  (ez=  0) 


eT(y,z)  =  ey(y)  e-jkzz  uz 
hT(i/,z)  =  hy(y)  e'jkzz  uy  +  h^(y)  e-jkzz  Uz 


kzYzhxCy/  z)  =  ®e 

Thus  (after  cancelling  e'jkz^  factor) 


lt  + 


Vx  Vx 

k2 


•  ^ 


zX  e  x(y,  z) 


kJxi^^y  +  h^^z> 


1 


mo 


der 


[ko£r-k^  jeyu^-jk^-^u 


If  is  chosen  to  be 


where 

then 


k| 


^  kzOJ^o 

1 2  ,2  ,2  ,2 

kx  +  ky  =  ko  er  -  kj 

tl  M 

^2  ^  ^y 

tt 

'  jkz 


‘y 


k|9y 
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APPENDIX  III  REGION  EIGENMODE  COUPLING 

In  this  appendix,  the  coupling  between  x-modes  in  adjacent  subregions  is  calculated  for  TE  and  TM-to-x 
and  z  modes.  The  derivation  is  essentially  similar  to  the  TE  and  TM-to-y  mode  coupling  derivation  in 
Section  2.3.2.2,  and  proceeds  using  the  same  approach.  The  x-mode  coupling  is  evaluated  by  considering 
the  continuity  of  their  constituent  y-mode  fields  which  were  derived  in  Appendix  II  for  TE  and  TM-to-x 
and  TE  and  TM-to-z  modes  and  Section  2.3.2.1  for  TE  and  TM-to-y  modes.  Note  that  complete  listings  of 
each  of  the  coupling  matrices  is  omitted.  Appendix  IV  contains  these  expressions  for  the  general  cases 
of  section  subregion  coupling  between  all  possible  combinations  of  TE  and  TM-to-x,  y,  or  z  x-modes. 


III.1  TE,  TM-TO-X  COUPLING 

IILl.l  Electric  Field  Matching 

To  evaluate  the  x-mode  coupling  at  subregion  interfaces,  the  components  of  the  x-modes' 
constituent  y-mode  electric  fields  tangential  to  the  interface  are  matched  across  the  interface. 

y-component 


Matching  the  y-component  of  the  electric  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 


limn  i.n  (OUn  imp  tp  ^  ami  »i  0}Ur\  ^  'J 


Using  the  orthogonality  condition  for  TE-to-x  y-modes  in  subregion  B 

dy  =  dij 


and  integrating  over  the  subregion  interface  yields  the  coupling  equation 


yu  N^'  ■  yu 

(OUn  S-i  'np  f  an  *  >P  -A  «mt  f  uti  4,  ni 

1  2'' A,  I  iB,(y)  U^^y'idy 


^  p  yi 


^  >  yi 


This  equation  can  be  expressed  in  terms  of  coupling  matrices  A,  B,  and  G  as 


Vb,  +  I  G„pVe^  =  X  +  X  KjVa, 


Note  that  A,  B,  and  G  developed  for  TE  and  TM-to-x  modes  are  equivalent  to  matrices  derived 
for  TE  and  TM-to-y  mode  coupling  at  subregion  interfaces. 
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z-component 

The  z-component  of  the  electric  field  is  matched  across  the  interface  between  subregions  A  and  B 
in  the  same  fashion  as  the  y-component  to  yield  additional  coupling  relations.  This  field 
matching  is  expressed  as 

yyTAf 

B  fj^p  B  „p  A  A  «mj  «j 

1  Vb,  ''b.  ''b, (^)  =  S  ''a. ''a, W - 7^  I  Va,  % (» 


Z  p 


,m 

;■ 


Using  the  orthogonality  relationship  for  TM-to-x  y-modes  in  subregion  B 

and  integrating  over  the  subregion  interface  yields  the  coupling  equation 


yu  ■  yu 

-T^  1  ^Byiy)  yBy(y)dy=  £  j 

p  i  yi 

y 

^  J  yi 

where  exactly  the  same  approach  is  used  as  for  Section  2.3.2. 2.  This  equation  can  be  written  in 
terms  of  coupling  matrices  F,  C,  and  D  as 


imi  rJL 


I  mn  -si,  »mp  i  mi  »«/ 

''b,  -  SF»Ab,  =  SC„,l^A,  -  2D„;''a 


III.1.2  Magnetic  Field  Matching 

To  obtain  the  remainder  of  the  coupling  expressions  needed  to  evaluate  x-mode  coupling  at 
subregion  interfaces,  the  x-modes'  constituent  y-mode  magnetic  field  components  which  are 
tangential  to  the  interface  are  matched  across  the  interface. 

y-component 

Matching  the  y-component  of  the  magnetic  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 


^  ,mn  <n  rrtcp.  ^  «mp  «p  ^  <mi  <i  (0£n  ^  ..mj  ,.y 


Z  p 


Using  TM-to-x  y-mode  orthogonality  in  region  A 

tyASy)yLiy)dy^5ij 
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and  integrating  over  the  subregion  interface  yields  the  coupling  equation 


tueo 


nI^  Ju 

^  ump  f  ’n  *  np 


M™  V 
^  ,n 


’a.  iv,(y)  V,(y)dy=  1  h,  j  ^yfv^  (y)dy 


yi 


yi 

,TE 


"TS-  E'e,  I  Vg  {>’)rfy 


^  J  y, 

This  equation  can  be  expressed  in  terms  of  coupling  matrices  <I>,  C,  and  A  as 


N 


TE  ^TE 

"”'p  .  ""y 


t  mu  ■  ■  urnp  ■  ■  ^  t  '■  -  i 

P  i  j 

z-component 

Matching  the  z-component  of  the  magnetic  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 


lu,  <.y)iA/y)+^2'A,  <.y)iA,(y)=  Ih,  W/g.O) 


m 


N 


TM 


mp  ,p 

y 


nF 


nF 


copQ  ^  j"y  .  ,  ,  , 

2  hyc  (y^^Byiy) 
H  j 


Using  the  TE-to-z  y-mode  orthogonality  in  subregion  A 

and  integrating  over  the  subregion  interface  yields  the  coupling  equation 


N™  .Vu 


nF  -.-y" 


nmn  frtLin  '^P  r  *  'P  r 

h S  J  'a  fw  'a  (W‘'y=  I V  I  wy>  Vy>^y 


yi 


Vi 

tTM 


ijri 

j 


yi 


This  equation  can  be  expressed  in  terms  of  the  coupling  matrices  T,  A,  and  P  as 


Ni 


Ni 


timn  imp  *  umi  in 

2  2  2  Pnj^B 


Nh 


imj 
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III.1.3  Complex  Power  Conservation 

As  for  the  TE  and  TM-to-y  mode  derivation  in  Section  23.2.3,  the  coupling  expressions  from  the 
field  matching  for  the  TE  and  TM-to-x  modes  are  written  in  matrix  form  as 


I 

1 

1 

-F 

G 

1 

1 

I 

I 

1 

1 

r 

1 

1 

I 

-D  V 


.tm 

A  JL^A. 


-A  I 


L  r  ;  I  jLi;j  L  p  1 jl^bJ 

The  power  flow  through  the  discontinuity  for  the  m*^^  z-mode  is  determined  from  the  fields  in 
both  region  A  and  region  B 

y.  ,  <<  ,yg 

=  1  E”,  X  H?,  <<)’  =  S  S  J  X  hf ^  dy 

Jl  n  p  yi 

yu  yu 

Pi  =  /  Efi^  X  Hfjdy  =  X  X  V^llf  I  ef.^  x  hf/rfy 


For  any  given  region,  ep  and  may  be  split  into  TE  and  TM  components  and  the  power 


- j  o- . o— -IB - iB 

integral  evaluated  as  follows 


Je^xh^  •Ux^/t=|  (eyUy+e^u^)x(hy  Uy+h^  uj  •u^dy 

yi  yr 

yu 

f  tfi  t  p  ifi  ip  jj, 

=  j  ey  (y)h^  (y)  (y)fiy  (y)  dy 

yi 

/rp'U.n  '«  'P  ±  in  <p  * 

(y)*,^+  J  (y)V,  (yfdy 


.mr  m  •' 
‘^z  yi 


,  tn  r  m  J  -y  -y 


,/n,m  •' 

‘^z  '^z  yi 
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Ju  *  yuf  *  *  ■ 

^  n/l  ftp  P  tin  ftW  tip  tip 

J  Cj'  xHt-  •Ux<i>'=  I  {CyWy+e^  Uz)x(/ty  Uy  +  ft^  u^)  •UT^dy 

\t.  yy.  ^  ^ 


»  tl«  H/;  5|,  ?  tin  Itp  ;(i 

=  J  Cy  (y)h^  (y)  dy-  I  (y)hy  (y)  dy 


y  2  ^ 

r  "ft  lip  if:  Oy  £«2  f  iin  lip  i|c 

j  ly  (y)/y  (y)  ^y+ — ^  J  (y)^j  (y)  dy 

J  •'  h-mrm  J  ^ 

yi  '^z  '^z  yi 


^5 

^np  ^  *  J 

'^z  y/ 


1*  It/t  tip  A 

J  Vy  (.yWy  (.y)  dy 


tirt  tp*  itn  tin  ip*  ip* 

J  Cy.  xhy.  •Ux^ij=J  (eyny+€^u^)x(hy  Uy+h^  u^)  •Uj^dy 

U.  L. 


yur  1 

^  tin  I  p  tk  '»n  I  p  )k 

J  “^z  ^y^S 


y^  in  lip 

J  e  J’  X  hj’ 


=  ^  j  /;  (y)ly  (yfdy  -  ^  J  V;’  ly)v'  (j)*  dy 
yi  ^  yi 

=  in  region  A  and  [g^  -  F  in  region  B 

ipn  L  Jpn 

y^  in  in  lip*  lip* 

'Uxdy=J  (eyUy  +  ezU^)x(hy  Uy  +  h^  Uj_)  •Ui^dy 

yr 

yur  -1 

=  J  ey  (y)h^  (y)  -e^  {y)hy  (y)  dy 

yi 

-^]'y  iy)'y  (>)*  <3’)<  w*  “y 


,m  •' 

'^z  yi 


F  — in  region  A  and  G  — in  region  B 
inn  inn 


Note  that  denotes  the  conjugate  transpose  of  the  superscripted  matrices.  From  the  completeness 
relationships  shown  below,  the  power  crossing  the  discontinuity  can  be  calculated  for  the  modes 
in  both  region  A  and  region  B. 
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So,  for  region  A 


tiF  „.2m?  y- 


r  +  1  ^  0)  Lir\  r  *  r  . 

G  G  =  ^  5r  J  J  ^ 

inn  ^  ,m,m  J  ^y  °y  J 


^  \  '^z  yi 
,.,2„2  . 


(O  tin  f  <n  *  .p 

= — ^  j  ig  h  ^y^^y 

,  m  ,  m  J  “v 


.m.m  •< 
'^z  ‘^z  yi 


rV]  =£ 


^  '^z  '^z  }'/  yi 


,m,/n  - 
\  '^z  yi 

<  ro2p2  y, 


/*  *n  *  '/2 


4.  1  69  Sf\  r  ♦  c  "P  * 

F^F  =  £  - ^  J  Vg  (y)  dy  J  Vg  (y'  )Vg  (y' )  dy' 

inp  ^  ,m,m  J  "y  ^y  J  ^y  °y 


imim  - 
«  '^Z  ^Z  >/ 

.,2„2  Ju  _ 


if"”  +  "P 


,m,m  •'  ‘"r 

'^z  '^z  yi 

Na^  „2„2  >-« 

=  y  ^  -%-f 

^  '^z  '^z  yi 


ru  yu 

r  »<?  »'W  ;jc  r  "P  + 

j  ‘'a,W%W  -iyj  (y)r,  (>  )  ‘‘y 


= 


9  9  V 

69  £r»  r  t*n  "P 

- ^A  (>’)  ^A 

.m.m  J  "r  '^y 

'^z  '^z  yi 


i+r+T  r+-®  V 


r-4>+  ,  i+$+o  V 


r+  i+r+r ,  r+-o  v 


r-$+  ,  i+o+o  V 
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From  the  relationship 


I 

-a. 

-1 

[i+on-i 

r 

1  I  - 

_-[i+r«)]-ir 

[i+ro]-i  _ 

'  ' 

-1“ 

1 

I 

[i+or]-i  I  «D[i+ra>]-i 


l_-r[i+or]-i  I  [i+r<i>]-i  J 

it  follows  that  [I+r«I>]'lr  =  r[I+<I>r]'l  and  [l+Or]"l$  =  ®[I+rO]‘l.  Using  these  new  relations 


i+r'*T 

r+-o 

I 

-o 

I  1  r+ 

r-o+ 

i+$t$ 

r 

I  _ 

I  I 

Hence  the  power  in  subregion  A  is  given  by 


■  - 

% 

t 

■c  1 

-V- 

--A*  1 

A 

"  ,m  - 

t 

C  +  B+G 

1  p+-C(b 

,fn 

_  Ar-A+ 

1  A+A+O 

Similarly,  for  region  B 


t 

I+G+G 

G+-F 

r 

Vr 

Bx 

..m 

Ir 

L 

_  G-F+ 

I+F+F_ 

,.m 

LvbJ 

"  ,m  - 
*Bx 

t 

I+G+G 

G+-F 

I 

-1 

C 

-d’ 

r  "1 

,fn 

_  G-F+ 

I+F+F_ 

G 

1  I  - 

B 

A  J 

,.m 

L^aJ 
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Using  the  relationship 


'  I 

-F 

-1 

[I+FG]-1 

[i+fg]-if 

_  G 

I 

_-[I+GF]-lG 

[I+GF]-1  _ 

'  I  1 

G+ 

-1“ 

1 

I 

[I+FG]-1  I  F[I+GF]-1 
-G[I+FG]-1  \  [I+GF]-1 


it  follows  that  [I+GF]-1G  G[I+FG]'l  and  [I+FGJ-^F  =  F[I+GF]-1.  Using  these  new  relations 


I+G+G 

G+-F 

"  I 

-F 

-1 

I 

G+ 

_  G-F+ 

I+F+F_ 

_  G 

I 

I 

Therefore,  the  power  in  subregion  B  is  given  by 


“  ,m  " 

T 

'  I  1 

G+ 

C 

1 

r  "1 

1 - 

' 

\ 

1 _ 

1 

I 

B 

A  _ 

..m 

L^aJ 

“  .m  “ 

t 

C  +  G+B 

1  G+A-D 

>  ' 
1 _ 

_  B-F+C 

A+F+D_ 

&  ^ 

!  > 
_ 1 

The  following  substitutions  will  ensure  that  for  any  number  of  x-modes  in  subregions 

A  and  B 

B  AT 
A->C+F 
P^A+G 
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As  in  Section  2.3.2.3,  these  relationships  must  be  satisfied  automatically  for  an  infinite  number 
of  modes  in  subregions  A  and  B,  since  both  the  field  matching  and  power  conservation  will  in 
principle  be  perfect.  This  can  be  proven  by  using  completeness  relations  for  the  y-mode  fields  as 
follows 

(i  Vi  Z  Vi 


Similarly 


VuVu  ^  „  ^'a 


2  vi  yi 


„ 


im  »  .f 


.m  J  Uy 

m 


B„p  as  required. 


[“4=  S 

1  « 


^  ViVi  ? 

=  -^  f  yl”  fwlV?  (v)dv 


2  yi 

jTM 

Z  fyWyiy' 


m  »  „p 


=  D„p  as  required. 

^B^y«  « 

A^G  =Xj;;fyj 

Jnp  " 

'  /7  1/. 


y;  2:  y; 

_  _Tr 


yuy« 

=»  f  f/; 

J  J  Aj 

2  yi  yi 


nF 

nfi  a-  fV  ^  *  >fq 

^y^  h  (y'^^h  ^y^  h  ^y''>^y<^y' 

y  y  y  y 


y^  ..n  .  ir 

I  >a(^>  h 


J  Ay 

z  yi 


(y>dy 


3„p  as  required. 
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Nl^Vu 


^  vi  ^  yi 


VuVu 

^  VlVi  ^ 


z  yi 

=  A„p  as  required. 


The  final  form  for  the  coupling  between  TE  and  TM-to-x  modes  in  adjacent  subregions  A  and  B 
which  ensures  conservation  of  complex  power  is  therefore 


I 

1 

1 

-F 

G 

1 

1 

I 

I 

1 

1 

-O 

r 

1 

1 

I 

C 

j  -C«I> 

AT 

A  , 

C+ 

-C+F 

A+G 

1  aV 

III.2  TE,TM-TO-Z  COUPLING 

IIL2.1  Electric  Field  Matching 

To  evaluate  the  x-mode  coupling  at  subregion  interfaces,  the  components  of  the  x-modes' 
constituent  y-mode  electric  fields  tangential  to  the  interface  are  matched  across  the  interface. 

y-component 

Matching  the  y-component  of  the  electric  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 


^  ..mn  itn  k  mUn  ^  <mp  ,p  ^  t,mi  ..i 

S  V'i),  I  v'b,  4,W=  I iA,(y) 

I^T  ^ 

^  m  P 


%  ; 

^  m  J 


Using  the  orthogonality  of  TE-to-x  y-modes  in  region  B 
rh  ,  "j  ,  *  , 


J‘rl  »V  * 

Qhy(y)hy(y)  dy  =  dij 
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and  integrating  over  the  subregion  interface  yields  the  coupling  equation 

yu  nI^  .yu 

vmn  k  rnitf\  -A  *nip  ?  nn  *  »P  ^  »rni  ^  .in  *  ui 

S  ''b.  I  isAy)  iB,(y)dy=  tv  A,  J  'A,(.y)‘‘y 

Kt 


p  yi 


N™  V 

X  vZ  I  hy(y)*  I  Ay  (y)dy 

I’m  j  yi 


This  equation  can  be  expressed  in  terms  of  the  coupling  matrices  G,  A,  and  B  as 

N™  Nf  .  N™ 

umn  ^  tmp  ^  tmi  ^  tm; 

Vb,  - 

P  i  J 

In  anticipation  of  the  further  derivation  to  ensure  power  conservation,  also  define 

^np=  \2  J  %  (3^)  Uy  (y)dy 


z-component 

Matching  the  z-component  of  the  electric  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 

Nl^  Nl^  .  . 

tvA,VA,(.y) 

n  i 

Using  orthogonality  between  TM-to-z  y-modes  in  subregion  B 

and  integrating  over  the  subregion  interface  yields  the  coupling  equation 

,  .i 

Vb,  =  X''a,  VA,<.y)‘‘y 

'  yi 

This  equation  can  be  expressed  in  terms  of  the  coupling  matrix  C  as 

Nl^ 

I  mn  I  mi 

Vb,  =  t^niVA, 


IIL22  Magnetic  Field  Matching 

To  obtain  the  remainder  of  the  coupling  expressions  needed  to  evaluate  x-mode  coupling  at 
subregion  interfaces,  the  x-modes'  constituent  y-mode  magnetic  field  components  which  are 
tangential  to  the  interface  are  matched  across  the  interface. 
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y-component 

Matching  the  y-component  of  the  magnetic  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 

,mn  ,n  k^^COEc)^  "’^P  "P  ^  ^  T/"-'  z'  \ 


‘m  P 


'm  J 


Using  orthogonality  of  TM-to-z  y-modes  in  subregion  A 


rh  'i  'j  * 


and  integrating  over  the  subregion  interface  yields  the  coupling  equation 


ajPE  V  N™  y 

^mn  k^COEn  ^  ?  in  :4c  "p  ^  r  /  \  j 


P  yi 


k^Caeo^iA  "mj^r  <n  *  "J 

^  1‘b,  I  ''a,w 

J  yi 


This  equation  can  be  expressed  in  terms  of  the  coupling  matrices  C,  and  A  as 


tmn  ^  nmp  ^  ^  imt  ^  vm 

^Ax  ^  ~  ^  ^in  ^B^ 


In  anticipation  of  the  further  derivation  to  ensure  power  conservation,  also  define 


F„p  =  J  v'b  (>)*  v'bI  iy)  dy 


z-component 

Matching  the  z-component  of  the  magnetic  field  at  the  interface  between  subregions  A  and  B 
yields  the  expression 


•Jx,  «mn  «n  ^  ><mi  ni 

S4,  /a,W=  S'b,  hAy) 


Using  orthogonality  of  TE-to-z  y-modes  in  subregion  A 


J'H  "J  jk 

QUy(y)Uy(y)  dy=Sij 


and  integrating  over  the  subregion  interface  yields  the  coupling  equation 


limn  — ^  nmi  /•  itn  ^ 

u,  =  I'b,  j  iAy(y)  iBJy)dy 
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This  equation  can  be  expressed  in  terms  of  the  coupling  matrix  A  as 


<mn  ^  II 


IIL2.3  Complex  Power  Conservation 

To  ensure  power  conservation  across  the  subregion  interface  regardless  of  the  number  of  x-modes  in 
subregions  A  or  B,  the  technique  described  in  Section  2.3.2.3  is  used.  By  constructing  column 
vectors  for  the  transverse  eigenmode  voltages  and  currents  for  the  m*^  z-mode,  with  the  TE  and 
TM  components  separated,  the  coupling  expressions  from  the  field  matching  can  be  written  in 
matrix  form 


I 

1 

1 

0 

-G 

1 

1 

I 

I 

1 

1 

4) 

0 

1 

1 

I 

The  power  flow  through  the  discontinuity  for  the  m*^^  z  mode  is  determined  from  the  fields  in 
both  region  A  and  region  B.  The  coupling  between  the  modal  fields  in  regions  A  and  B  must  ensure 
that  the  power  coupled  from  modes  in  region  A  to  modes  in  region  B  is  the  same  as  the  power 
coupled  in  the  opposite  direction. 

Pi  =  jEf,  X Hf,  I  'h  •'r.  ‘‘y 

yi  P  yi 

yu  yu 

Pi  =  J  E?,  X  HfJ'dy  =  X  S  j  X  b^Jdy 


For  any  given  region,  and  may  be  split  into  TE  and  TM  components  and  the  power 

integral  evaluated  as  follows 

r  _ 

?  in  ip  f  <n  I p  np 

JeyXhj.  •Uxd>’=  J 

yi  yr 

=  I -e^  iy)hy{yfdy 


=  J  y'ym'y  iyfdy 
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n  n* 

ll?l  lip 


>'«r 

I*  iin  lip  ■  /• 

I  67’  X  hj  •  Ux  =  J 


yi 


ttfi  ttfi  fip 

(.^y  Uy  +  Uj )  ^ 


J/ 

>’«r 


•  Uj,dy 


=  I 


iiM  i»P  ^ 

ey  (y)K  (y) 


dy 


yi 


y« 


/•  iiW  lip  4. 

=  -\ly  (y)Iy  (y)  dy 


yi 


=  d, 


np 


y«  yu 

e  '«  <'P  r 

Cy  Xhj’  •Ux^y=  I 


yi 


yi 

yu 


tn  in 


(eyUy+eyU^)X(hy  Uy+k^  ) 


•  u^dy 


»•  ttl  lip  4,  in  ''  y 

=  J  -^y  (y)/!^  (y)  (y)/iy  (y) 


y« 

2’ 

"r  yi 


n 

tl/l  I  P 


p  tl/I  I  P 

I  67’  X  hj’  •  Ux  4y  =  0 


yi 

Using  the  previously  defined  TE  and  TM-to-z  z:-mode  coupling  matrices,  the  power  in  region  A  is 
given  by 


= 


,.m 


L^BxJ 


,in 

Ir 

Ox 

,fn 

L^xJ 


,>m 

^B*J 


I  1 

0 

o+-r  1 

I  J_ 

c  [ 

0 

1 

A  J  _ 

c 

!  0- 

A+-A<D+ 

1  A. 

c 

0 

A+-Ar  1 

A  J 

,fn 


$+  ,  I 


O+'F  I  I 


,pn 


<D+-r 

.m  ■ 

L^AxJ 


I  JL^a 


.m 


...m 

x- 


V, 
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Similarly,  for  region  B 


Note  that  denotes  the  conjugate  transpose  of  the  superscripted  matrices.  From  the  completeness 
relationships  shown  below,  the  power  crossing  the  discontinuity  can  be  calculated  for  the  modes 

in  both  region  A  and  region  B.  The  following  substitutions  will  ensure  that  for  any 

number  of  x-modes  in  subregions  A  and  B 

A  -^C+F 
B  ^AF 

As  in  Section  2.3.2.3,  these  relationships  must  be  satisfied  automatically  for  an  infinite  number 
of  modes  in  subregions  A  and  B,  since  both  the  field  matching  and  power  conservation  will  in 
principle  be  perfect.  This  can  be  proven  by  using  completeness  relations  for  the  y-mode  fields  as 
follows 


^  as  required. 

inp 
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Similarly 


N'^yu  n 

„  I*  vn  *  iio  K,  (OUo  t  'P  "9  it 

lAr4=  1 1  /b,W 


<l  Ji 


J/  yi  9 


k. 


yuyu 


k^j 

,m. 


'm  J/ 


=  ^72— J  '«  (>■>  'a  (>’>‘'2' 


'm  y/ 


=  Bfip  as  required. 

The  final  form  for  the  coupling  which  guarantees  conservation  of  complex  power  regardless  of 
the  number  of  x-modes  in  subregions  A  and  B  is  therefore 


I 

1 

1 

0 

-G 

1 

1 

I 

I 

1 

1 

0 

r 

1 

1 

I 

,m 

,.m 


..m 

L*A;,J 


-AF 

C+ 


0 


0 


A 

C+F 


A+ 


,m 

l^a,. 

.m  ■ 

,.m 

L^Bj:J 
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APPENDIX  IV  SUBREGION  COUPLING  SUMMARY 

The  coupling  between  adjacent  subregions  has  been  determined  using  a  field  matching  technique  in 
Section  2.3.2.2  for  TE  and  TM-to-y  y-modes  and  Appendix  III  for  TE  and  TM-to-x  and  TE  and  TM-to-z 
y-modes.  The  coupling  expressions  subsequently  obtained  ensure  that  the  complex  power  crossing  the 
subregion  interface  was  conserved,  while  still  satisfying  the  field  matching  criteria  when  an  infinite 
number  of  modes  are  included.  The  coupling  expressions  for  two  subregions  where  the  x-mode  fields  are 
expressed  in  terms  of  mode  functions  which  may  be  TE  and  TM-to-x,  y,  or  z  are  listed  in  this  section. 
The  full  derivation  is  not  considered  worth  repeating,  and  is  almost  identical  to  the  approach 
described  in  Section  2.3.2.2  and  Appendix  III.  The  full  coupling  expressions  for  arbitrary  y-mode 
reference  directions  allows  the  use  of  different  combinations  of  TE  and  TM-to-x,  y,  and  z  y-modes  in  the 
same  two-dimensional  cross-section. 


IV.l  TE,  TM-TO-Y  MODES  IN  A,  TE,  TM-TO-Y  MODES  IN  B 
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IV.2  TE,  TM-TO-X  MODES  IN  A,  TE,  TM-TO-X  MODES  IN  B 


c 

1  -C$ 

AF 

1  ^ 

C+ 

1  -C+F 

A+G 

1 

/*  tm  ^ 

Kp  =  \hy(y)  iAy(y)dy 

yi 

C^  = 


^'.p=^!wy^’''liy)dy 


’»pp=^iWy'>'''A,^y^<‘y 


^pp  =  ^!'Bl(y)'‘fB^(y)dy 
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IV.3  TE,  TM-TO-Z  MODES  IN  A,  TE,  TM-TO-Z  MODES  IN  B 
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IV.4  TE,  TM-TO-Y  MODES  IN  A,  TE,  TM-TO-X  MODES  IN  B 


G 

I 


L^b,. 
^ ' 


0 

C+ 


A+G 


-CO 


A 

-C+F 


A+ 


L^A,. 

.m' 

lVj 


^np  -  J 


yuig  (yfiAiy) 

'  - ^dy 


yi 

yu 


£A(y) 


yi 

<>■)*<  W”* 


Z  yi 


k’^mo  'f  *  «p 

I  y(y>  \W‘'y 


^'^pa  yi 

z  y, 


np  -  ,2 

yi 


£A<<y) 
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IV.5.  TE,  TM-TO-X  MODES  IN  A,  TE,  TM-TO-YMODES  IN  B 
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IV.6  TE,  TM-TO-X  MODES  IN  A,  TE,  TM-TO-Z  MODES  IN  B 


-G 

I 


I 

-<[) 


I 


■ 

L^B,. 


.tin 


AE 

C+ 


0 


-CO 


A 

C+F 


yu 


A+ 


<•  n7I  ^  »P 

Kp=\hy(.y)  ^Ay(y)dy 


.m 


,fn 


L^J 


yi 

yu 


C„p^jVB(.yfvUy)dy 


yi 


/•  3lj  tiO 

j  ^B^y)  yByiy)dy 


yi 


yi 


,  k^COjilQ^r  JP  ,  .  , 

^np  =  ^72^J  hAy)  %{y)dy 


yi 


yi 
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IV.7  TE,  TM-TO-Z  MODES  IN  A,  TE,  TM-TO-X  MODES  IN  B 


.m 

,4n 

l^a,. 

fa  ^ 
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-F 

r  .in" 
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1 
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L^b,. 

_-Ar 
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I 
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4) 

r 

c+ 

-C+F 
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I 

L^aJ 

_A+G 

A+ 

JL^b^J 


/•  »i«  ^  vp 

Kp  =  J  hy(y)  Uy(.y)dy 

yi 

^np  =  \yBy{yf^'lyiy)dy 
yi 

Z  yj 

■^n,-^hl(yryZ(y^‘^y 


yi 


yu 


G.p=^lhAy)'iliy)‘^ 


yi 

yu 


^np=  \2  J  ^Ay  (y)  I  Ay  (y)  dy 

Kt  ^  ^ 

y, 
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IV.9  TE,  TM-TO-Z  MODES  IN  A,  TE,  TM-TO-YMODES  IN  B 
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-F 
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..m 

l^b,. 
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-Ar 
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I  %(y'>  iA/y)‘iy 
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C„p=  jvB^(yfV^^(y)dy 

yi 

Wvi'wrfy 


‘^pb  yi 

m _ _  yu 
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APPENDIX  V  ANALYTIC  ZERO,  POLE  LOCATION  TECHNIQUE 

From  the  moments  Sn  derived  from  Cauchy's  Theorem  as  discussed  in  Section  2.3.27,  the  position  of 
zeroes  and  poles  of  a  function  within  a  given  contour  can  be  determined  analytically  in  the  following 
cases. 


So  =  0 

Can  handle  no  zeroes  or  poles,  or  one  zero  and  one  pole.  If  there  are  no  zeroes  or  poles,  then  Si  also 
equals  zero.  For  one  zero  and  one  pole 


5i=zi-/?i 

*^2  =  -  Pi  =  (zi  -Pl){zi  +Pl)  =  Si{zi+p^ ) 


So=  1 


Pi 


1 

2 

J_ 

2 


Can  handle  one  zero,  or  two  zeroes  and  one  pole.  For  one  zero,  zi  =  Si.  For  two  zeroes  and  one  pole 


So 


and 


5i  =zi  +Z2-P1 
^2  =  Zi  +Z2-  Pi 
S'i=zi  +Z2-P1 


Z1+Z2 


2 

iS2-Si 


Pi  =  C- Si 


Hence,  zi  and  Z2  may  be  separated  as 


21  = 


22 


C  +  2Sl  -  4CSi  +  IS2  j 

+1S'1  -ACS1  +  IS2 


101 


DSTO-RR-0027 


So  =  -l 

Can  handle  one  zero,  or  two  poles  and  one  zero.  For  one  zero,  pi  =  -  Si-  For  two  zeroes  and  one  pole 


So 


and 


5i=Zi -P2~P\ 
c  .,2  2  1 

■^2-^1  -P\  ~P2 

c  ..3  3  „3 

■>3  =^1  ~Pl  ~P2 


Pl+Pl  = 


3  ^2  +  Si 


=  C 


Zi  =  C  +  5] 


Hence  pi  and  p2  may  be  separated  as 

Pl=-  C  +  ^Jc^+2S^+4CSl-2S2 

2  _ 

P2=|  C-^Jc'^+2Sl  +  4CSi-2S2 

So  =  2 

Can  handle  two  zeroes 


where  zi  and  Z2  are  given  by 


So  =  -2 

Can  handle  two  poles 


where  pi  and  p2  are  given  by 


5i=Zi  +Z2 
^2  -  ^2 


Si=Pi  +P2 

C  rt^  -L 

^2- Pi  +  Pi 
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APPENDIX  VI  CROSS  SECTION  FIELD  EXPANSIONS 

For  evaluating  the  z-mode  fields  in  each  section  for  output  and  for  calculating  the  coupling  between 
z-modes  in  adjacent  sections,  the  z-mode  fields  are  expressed  in  terms  of  the  A:-mode  and  y-mode 
equivalent  network  voltages  and  currents.  These  voltages  and  currents  are  defined  in  Section  2.3.1  for 
y-modes  and  Section  2.3.2  for  x-modes.  The  derivation  of  the  x-components  of  E  and  H  from  the  y  and 
z-components  is  described  in  Section  2.3.2.9.  The  y  and  z-components  are  written  directly  from  the 
definitions  of  the  y-mode  and  x-mode  field  in  Sections  2.3.1  and  2.3.2.  The  power  transmitted  in  the 
z-direction  by  the  z-modes  in  each  two-dimensional  section  is  derived  from  the  field  expansions  as 
discussed  in  Section  2.3.2.9. 


VI.l  TE,TM-TO-YMODES 

jTM  „n 


N, 


n 
nTE 

•‘’V  V  7  m 

^  ^  imn  in  ^  »mn  un 

ET  MVy  (,)-  2  (x)V^  (.y) 


n  Kn 


TE 


'  ,mn  ,n 
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NT'S  j^TM 

Hz  =  E  ^"-^4  Miy  ()’)+  E  4  W4  <y) 


n  Kn 


Therefore 


nTE  j^TM 

=  I  (^)Vy  (y)-  X  4  MVy  (y) 


'un 


and 


NEE  f^TM 

\  1  TYl 

^  yJL  tmn  *n  ^  k  Ci)Pf\ 

Hx  =  t.''x  UK,  (y)- UKjW 


'•Mn 
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The  power  flow  in  the  z-mode  through  the  two-dimensional  section  is  given  by 

„  w— 1  ■J-,  r  ""in  <mp  *  /•  nn  ip  * 

^Z  =  S  S  j  4  (^)  dxjVy  iy)Vy  (y)  dy 

n  P  X[ 


Therefore 
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^  ^  X  «mn  imp  ^  ly  (y)/y  (j) 

i  MV,  Mdxl^ - — dy 


n  P  xi 
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vt  t  j/rw/rwAj<w<w> 

«  P  ^un  xi  yi 

jyTM  jyTM  *  X  V  x""  r'P  * 

(OEn  r  ""1”  *  C  (y)^V  ^y^ 


n  P  "^up  xi 


VI.2  TE,TM-TO-X  MODES 
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The  power  flow  in  the  z-mode  through  the  two-dimensional  section  is  given  by 
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VI.3  TE,TM-TO-Z  MODES 


■*’v  »  Wl 

^  ^  .tmn  tin  ^  k  fDLLn  »'* 

=  S''*  (J’)-  S  W 


Af, 


rt 


^  '  I  mn  in 

E?  =  S  U)V',  W 


AT™ 

^’y  T  /w 

_  z  •  mw  t  n  _  Z  /,»'*  y 


^  i  imn  tn  ^  k^Cl)Pr\ 

ff”  =-  S  W-  S  ^3-“-^*  (>’) 


A^, 


n 

TE 


^  _,  timn  fin 

W/yW 


eT=-U-i,  (*»';,  W 


n  ^m 


COEqE^. 


^  ,mn  in 

=  Z  W4  (y)-  Z  W4  (y) 


fy/io 


105 


DSTO-RR-0027 


The  power  flow  in  the  z-mode  through  the  two-dimensional  section  is  given  by 
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APPENDIX  VII  SECTION  COUPLING 

To  evaluate  the  coupling  between  2-modes  in  adjacent  sections  A  and  B,  the  components  of  the  electric 
and  magnetic  fields  tangential  to  the  section  interface  were  matched  over  the  interface  The  technique 
used  for  the  derivations  in  this  appendix  is  identical  to  the  approach  used  in  Section  2.3.3.1  for  TE  and 
TM-to-y  y-modes  in  both  section  A  and  section  B. 

Calculation  of  the  section  z-mode  coupling  allowed  the  individual  sections  to  be  cascaded  together  to 
form  the  complete  element  being  analysed  as  discussed  in  Sections  2.3.3.2  and  2.3.3.3. 
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In  addition,  multiplying  both  sides  of  Equation  VII.l  by 


and  integrating  over  the  aperture  yields 


where 


^mn  ^  ^ 
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k 


(VII.2) 


Both  sides  of  Equation  VII.2  were  multiplied  by 

ivF  ^^b-K'I 

1^= - ^ ^ 

5  (O^O 


rs  »'r.y  4.  nj  ^ 

-^VB^^fh^^y) 


Integrating  over  the  aperture  in  the  x  and  y  directions  and  using  orthonormality  of  the  y-modes 
yields 


HV«  =JV, 
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where 
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VII.1.2  Magnetic  Field  Matching 


To  evaluate  the  z  -mode  coupling  between  adjacent  sections  A  and  B,  where  the  y-modes  are  TE 
and  TM-to-A:  in  both  sections,  the  x  and  y-components  of  the  magnetic  field  tangential  to  the 
section  interface  are  matched.  This  section  of  Appendix  VII  describes  this  process,  which  closely 
follows  the  TE  and  TM-to-y  mode  analysis  in  Section  2.3.3. 1. 


:c-component 
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Both  sides  of  Equation  VII.3  were  multiplied  by 


<,pq  «q 

(x)lByiy) 


(VII.3) 
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TE 


’Zv.ixfhly)' 


Integrating  over  the  aperture  in  the  x  and  y  directions  and  using  orthonormality  of  the  y-modes 
yields 


KI.4.=L*«, 
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where 
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In  addition,  multiplying  both  sides  of  Equation  VII.3  by 
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C0{1Q 
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and  integrating  over  the  aperture  yields 
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(VII.4) 


Both  sides  of  Equation  VII.4  were  multiplied  by 


A^™  ko^A-^xA. 


1.2  =' 
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Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the  y- 
modes  yields 


cia,  =m«. 
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where 
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VILI.3  Complex  Power  Conservation 

As  for  TE  and  TM-to-y  modes,  using  the  field  relations  in  Section  2.3.2.8,  the  power  flow  through 
any  given  section  is 
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In  terms  of  the  previously  derived  coupling  matrices 


II 

t 

d  +F^  -bK^] 

'''aJ 

= 

t 

and 


^[A  +  H-S][Vb_' 

^[b+j-t][vaJ 

So  for  conservation  of  complex  power,  require 

dUn^-l^=b+j-q 

This  is  true  by  inspection.  Therefore,  to  ensure  complex  power  conservation  regardless  of  the 
number  of  z-modes  used  in  each  section,  the  following  expressions  for  the  current  and  voltage 
coupling  will  be  used 

[Ci]  =  [C  +  F  +  Kr^  [D  +  G  +  L] 

[Cv]  =  [A  +  H-Sr^  [B  +  J-T] 

=  [A  +  H-Sr^  [dUg'^'+L’^' 


VII.2  TE,TM-TO-Z  COUPLING 

VII.2,1  Electric  Field  Matching 

To  evaluate  the  z  -mode  coupling  between  adjacent  sections  A  and  B,  where  the  y-modes  are  TE 
and  TM-to-z  in  both  sections,  the  x  and  y-components  of  the  electric  field  tangential  to  the 
section  interface  are  matched.  This  section  of  Appendix  VII  describes  this  process,  which  closely 
follows  the  TE  and  TM-to-y  mode  analysis  in  Section  2.3.3.I. 


112 


DSTO-RR-0027 


Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the  y- 
modes  yields 
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In  addition,  multiplying  both  sides  of  Equation  VIL5  by 
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and  integrating  over  the  aperture  yields 
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i/-component 
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Both  sides  of  Equation  VII.6  were  multiplied  by 
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Integrating  over  the  aperture  in  the  a:  and  y  directions  and  applying  orthonormality  of  the  y- 
modes  yields 
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In  addition,  multiplying  both  sides  of  Equation  VII.6  by 


nI^ 


and  integrating  over  the  aperture  yields 
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where 
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VII.2.2  Magnetic  Field  Matching 

To  evaluate  the  z  -mode  coupling  between  adjacent  sections  A  and  B,  where  the  y-modes  are  TE 
and  TM-to-z  in  both  sections,  the  x  and  y-components  of  the  magnetic  field  tangential  to  the 
section  interface  are  matched.  This  section  of  Appendix  VII  describes  this  process,  which  closely 
follows  the  TE  and  TM-to-y  mode  analysis  in  Section  2.3.3.I. 
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J  (VII.7) 


Both  sides  of  Equation  VII.7  were  multiplied  by 

Nf 

Svw  AW 


Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the 
y-modes  yields 
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where 
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In  addition,  multiplying  both  sides  of  Equation  VII.7  by 
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Both  sides  of  Equation  VII.8  were  multiplied  by 
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Integrating  over  the  aperture  in  the  x  and  y  directions  and  applying  orthonormality  of  the 
y-modes  yields 
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In  addition,  multiplying  both  sides  of  Equation  VII.8  by 
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VIL2.3  Complex  Power  Conservation 

As  for  TE  and  TM-to-y  modes,  using  the  field  relations  in  Section  2.3.2.8,  the  power  flow  through 
any  given  section  is 
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S  S  j^;c  W  ^\^y  ^y 


P  q  xi 
ritm  ‘'y 


yi 

yu 


''v  ‘’v  M  ■'U 

K  k  f  '"9  ""'P  *  e  '1  "P  * 

+  ^  I  S  S  ^\^y  ^y^^y  ^y'>  ^y 


P  q  Xi 

k  k  ^  ^  p  ,nq  »mp  *  p  op  ^ 

X  S  w  ^l^y  ^y^^y  ^y^  ^y 


yi 

yu 


'n  P  q  Xi 

n  i'ly  i'^y  X^ 

r  * 


yi 

yu 


+^Z  Z  I',  W',  w*&jv,wv,  w> 


k^k!}cOLLo£r  ^  ^  r  <nq  <mp  *  .9  .p  * 

+  ,  V  Z  Zl''.  W''.  W’*|',&U,W  rfv 


p  q  Xi 

N™N™x, 


yi 


VTm  P  9  Xi 

n  *  mTE  j^TE 

'.Xim  P'y 


yu 

!■ 

y/ 


t2 


knk^  COEn  f  "n^  "Utp  ^  «q  «p  * 

+  w,  y  z  z  j'«  W4  w*jv,  wv,  w'ij 

yi 


k  f  "”9  ""'P  *  f  ”P  * 

+  ^  X  X  l^;c  W  w^y  ^ 


>/ 


p  q  xi 

In  terms  of  the  previously  derived  coupling  matrices 

^4=[lAj[c^-F^  +  K^  +  Mt][v^^^ 


=  I 


B, 


D^-G^L'*'  +  N^||V 


][' 


‘2  J 


and 


f'?=K][A  +  H-P  +  S][v,J 
=  [l.,f[B  +  J-Q  +  T][v^_' 

So  for  conservation  of  complex  power,  require 

D^-G^  +  L^  +  N^=B  +  J-Q  +  T 
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This  is  true  by  inspection.  Therefore,  to  ensure  complex  power  conservation  regardless  of  the 
number  of  z-modes  used  in  each  section,  the  following  expressions  for  the  current  and  voltage 
coupling  will  be  used 

[Ci]  =  [C-F  +  K  +  Mr^[D-G  +  L  +  N] 

[Cv]  =  [A  +  H~P  +  Sr^  [B  +  J-Q  +  T] 

=  [A  +  H-P  +  Sr*  [d'^-GUlUn^' 
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APPENDIX  VIII  SECTION  COUPLING  SUMMARY 

The  coupling  between  adjacent  sections  with  the  same  choice  of  eigenmode  (TE,  TM-to-x,  y,  or  z)  has 
been  determined  using  a  field  matching  technique  in  Section  2.3.3.1  for  TE  and  TM-to-y  y-modes  in  both 
sections  in  the  element,  and  Appendix  VII  for  TE  and  TM-to-z  or  z  y-modes  in  both  sections  in  the 
element.  The  coupling  expressions  are  formulated  in  such  a  manner  as  to  ensure  that  the  complex  power 
crossing  the  two-dimensional  discontinuity  is  conserved,  regardless  of  the  number  of  z-modes  that  are 
included.  The  coupling  expressions  for  two  sections  A  and  B  where  the  fields  are  expressed  in  terms  of  y- 
modes  with  arbitrary  reference  directions,  i.e.,  TE,  TM-to-z,  y,  or  z,  are  listed  in  this  section.  The  full 
derivations  were  not  considered  worthy  of  repetition  and  follow  Section  2.3.3.1  and  Appendix  VII. 


VIII.l  TE,  TM-TO-Y  MODES  IN  SECTION  A,  TE,  TM-TO-Y  MODES  IN  SECTION  B 

[Ci]  =  [C-K  +  Mr^  [D-L  +  N] 

[Cv]  =  [A  +  H -  P]"^  [d^  -L^  +  N^ 


k 


— S —  J  M  dx 


Amn=l 


‘^Pb 

p  unq  I  mp  ^  P  '  P  * 

9  Xi  yi 


Hntn=  £ 


B  P  Q)£q  np  . 


f  «np  «mp  * 

J  Vb^  (x)Vb^  (x)  dx 


P  '^Upg  Xi 


.n,  ,n.p  *  ^viByiy)iBy(<yf 
Pmn=^  £  l^B^  ix)VB,  W  dx  j  ,  \ - dy 


p  q  xi 


^  ^  COjlQ  .nq  <mp  *  <q  ,p  * 

P>mn  =  X  T  dx  fVBy(y)VAy(y)  dy 

P  q  Xi  yi 

,nq  .mp  ,  ,q  .p  , 

^mn  =  £  £  j^Bjx)lB^  (x)  dx  j  VB^(y)VB^(y)  dy 
P  q  Xi  yi 

,nq  .mp  ,  il(y)iAy(y)\ 


£B(y) 


^  XI 

^md  'TkM 


M™  « 
/ya  i.n 


P  k^^(0£Q  u 


r  k 


^Aiy) 


f  .nr  .mp  *  (Y)^A,  (Y) 

I  Vb  ix)VA  (x)  dx  [  ^ - 

}  B,y  J  A,  K  J  J 


Ur  g  Z; 
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^  k,  COUq  e  ‘np  <mp  * 

Cmn  =  S  ,2-"-  I  (^)^A,  (^)  ^ 


^  '^“Pa 


^mn=  S 


<•  n/zo  vrnu 

J  ^A;,  (^)^A^  (^)  ^ 


L  “Pa 

....  *  .,  ...  * 

Mmn=t.t.  (a)  ^J^A,W^A/y)  ^X 

P  R  xi  yi 


VIII.2  TE,  TM-TO-X  MODES  IN  SECTION  A,  TE,  TM-TO-X  MODES  IN  SECTION  B 

[Cl  ]  =  [C  +  F  +  K]“^  [D  +  G  +  L] 

[Cy]  =  [A  +  H  -  Sr^  [d^  +  +  L"^ 


‘i. 


-'-AB  ,np  ,mp  . 

^  ^  k"‘<oe  e 


r  ^  9  * 

O  HA 

®  ?  ..n.  ..m.  ^ 

r.  ^ - S - -  J  ^B^  ^ 

2^  r  _2  *■ 

+  X  ^ - * - ~i^B^(^WB^  W  ^  ^>’ 


pm  *pn  J 

‘^Zb  '^Zb  xi 


n¥nW  knEB-k'l  Vu 

^».»=  S  S  „  W**Jv*^WVs^(y)  dy 

p  q  ^zb^zb  ^b  xi 


Dmn=  I 


^B 


kn^A  ~  kxA 

f  .  JT/g(<W4v*(y)v;(y)Vy 

k^  coeqEa  *' 


"  fc'”  Ea 
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<^mn  =  S  X  - - W  dx\lB  ij) I aAj)  dj 

■“  ■“  .M  itn  J  X  X  J  y  y 


“  “  un  im^  J  ^  ^  J  >"  y 

P  <1  '^Zb'^Za  xi  yi 

^^b~Kr 

^  D  ^  D  f  »nq  »mp  ^  ^  »q  »p  ^ 

^'«n=  I  2^1^: - ^jVBjxW^^  (xfdxjle^i^yM^^iyfdy 

p  q 

T  2  '2  *' 

^0  ~  ^xA  ^u 

^  ^mp  f  I W  imp  ^ 

^TM  - -i‘xJ’‘'>’A,  w  * 

Ja  t”  moEA  ,, 

^mn  “  ^ 

-£  - „  Jv  MIa^  M  dxjv^^iyW^^iy)  dy 

«  ^za  ^za^^  ^i  yi 

N™  Na^  knEA  -  k'l  Xu  yu 

^mn  =  X  £  ^ - i - -  (2C)V^  (x)  J /.  (y)/.  ( j)  dy 

P  <1  '^Za'^Za  xi  yi 

knEA  -  K  a  -'u 

^mn  X  ,n  J  M  dx 

p 


VIII.3  TE,  TM-TO-Z  MODES  IN  SECTION  A,  TE,  TM-TO-Z  MODES  IN  SECTION  B 

r/~<  l  —  f/''  1?  1  1  iv/t1~1  fTk  i''  1  T  1  vrl 


[Ci]  =  [C-F  +  K  +  Mp  [D-G  +  L  +  N] 
[Cv]  =  [A  +  H-P  +  Sr^  [d‘^-G^+l’'+N^ 


Nm  (DEqEb 


^7d  r  •'^p  '"'p  * 

]  ^B^  (x)Ib  (x)  dx 

J  X  ^x 


^tnn  ^  ,,Ti7 

I  2  ^ H 

P  kA  aA  A  Iinq  imp 


■n  ^  r  *^nq  ^mp  *  r  'P  * 

^  £  J/fi,  W/fi,  (xfdxjVB  iyWB^iyfdy 


\  Q  Xi 


«mn='L 


nI^  CO^O 


f>  unp  ump  ^ 

I  ix)VB^  (x)  dx 


7  k^*k^  ..„ 

~~T~^  5L  \'^B^ix)VB^  {x)*dx\lB  {y)lB  iyfdy 
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^mn  ~ 

P 


S  1  (x)*dx\ i"b^ iy)lBy {yfdy 


^Tb„  ^  xi 
Nr^  lrl-h.fl 


yi 


r  ^Tb„  Xi 


»^jV;:(x)v'«:wcfa|/;(,)4;(y)V, 


yi 


Nf 


^mn  ^ 
P 


j^n  i^m*  N™Xu 


S  jlBrU')h7(xfdxlvl(.y)vl(yfdy 


^Tb^  ^B  q  X[ 

khz* (Oe^ 


+■ 


kpB^  kpB^  r  Xi 


yu 

S  „nr  «mp  ^  f  nr  «p  ^ 

X  jh,(x)iB,  (x)  dxjvB^iyWB^iy)  dy 


yi 


N™ 


k-^mn  X 

P 


m  *  AjP^  X 
itn  A'g  x^ 


0)£q£a 


^  ,nq  .mp  *  ,q  ,p  * 

£  U)'&Jv'fi/y)V'A,W  dy 


q  xi 


yi 


+■ 


j.m^,n  N‘b‘^x^  _  __  Ju 

'^Za  ‘^Zb  ^ 


kpB^^A  r 


£  W/a^  (x)*d!xJyB^(y)V^^(y)*<<y 


yi 


N™ 


Gmn  -  X 

P 


'm*.n 


rm  jM  lyh  u 

1^7 .  'A  r  "'*<?  •  '"P  jk  r  * 

£  J‘'b,  W&  j/B^(y)/A,{y)*<<y 
*rA  ’  *'  « 

knCOjlQpB  :ic  ’''  'P  ♦ 

S  1^b,WFa,  U)**:J/B/y)/A/y)*<<y 


*rA  *rB  ''  ’I 


yi 


Nl^ 


k'mn  X 

P 


,n 


CUjUo 


Vr  -^M  3'h 

«  f  ""9  "ntp  *  /•  .t?  i.p  * 

X  \^B^  w  ^  dy 


q  xi 


yi 


T  N™X 
'2^ 


k^Pn  ^  r  nmp  ^  ^  tr  up 

"tV  X 

r  Xi  V; 


w„„  =  £ 

P 


1 2  Xu 

ko  ^  f;"9.  X* 


J 

y/ 

3'« 


r  '^9  :1c  r  '9  * 

£  W*JVi,^(y)VA/y)*rfy 


■^m 


y/ 


Zr2^W  /x-ir*  Xu 

kQk^  0}£q  "  (.  .,nr  -/«P  *  r  -r  ,.p 

+ - ! -  X  J  ^B/y)^A/>') 

yi 


yu 


krl  kri  r  x, 


'^TBu 
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VIII.4 


AT™ 


Cmn=  I 

P 


'^Za  f  -  P 


(OeQEA 


r  f  ftp  tmp  * 

J  (x)Ia^  (x)  dx 


+ 


].m*<n  nI^^u  yu 

I  JC^C M’dxjv7(y)vl(yfdy 


^r™ 


^mn  “  ^ 
P 


Na^^u 


y/ 


yu 


£  (xWa,  W  dxjlA^(y)lAy(y)  dy 


?  XI 


)'/ 


kQCOIlQEA^A  s*  ,  ■f/''  .  .j'P  . 

— — s - X  i^A.M^A,  w  dxjlA^(y)lA^iy)  dy 

^tI  kA  r  XI 


yu 


yi 


Nf 


Kmn=X 

P 


ift 

,  P  unp  itmp  a. 

— ^  J  (^)  dx 

CO/Iq  j  ^  ^ 


xi 


X  J  (^)^a7  (^)*  dx  j  I'Z  (yvly  (y)* dy 

^A„  q  Xi  V, 


ATf 

^mn  ~  Xi 

P 


jt2  ,nq  „mp  *  ,q  ,.p  * 

^  X  1^4  W'a,  W  aJva  (y)^^  (y) > 


‘^4  «  '=( 


^O^Za^^O  ^  f  /'"'■/  Nr"'”'’/  N*j  ^Ci/"'  I'  ^*J 

■  ^  t  ^  Z  J  W^A^  ^  J  ^A/>')''^A/y) 

^^A^  ^Ta„  r  xi  yi 


TE,  TM-TO- Y  MODES  IN  SECTION  A,  TE,  TM-TO-X  MODES  IN  SECTION  B 

1-1 


[Ci]  =  [C  +  M-Kp  [D-L  +  N] 
[Cv]  =  [A  +  H-Sr^ 


D^-L'*'  +N'^ 


Nr 


Amn  =  X- 

P 


2  '  ^ 

h^B-Ks 


—  J  {xf  dx 
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2  "2 
ko^B-kxB 
L - ^ - Jfy 

COHO  i  ^  ^ 

1/  _  V  2b 

^/nn  - 


Lin  *Kn  J 

Zb  % 


/•  tna  t.mp  ^  C  <,p  ^ 


=  £  S 


2  '  ^ 

koEB-^xB^q  , 


f  ina  nmp  *  t  <Q  "P  * 

I  (^)^B^  (^)  ^  J  '^B^  (>')^B^  W 


t"  t"*  c 

P  q  '^Zb'^Zb  Xi 


<  khT* 


^mn-  1,  7-^  . 


p  k^  B)]Ho  '”'■  '"’B  *  ■'■  'B  * 

_i_  V  _S4 _ ^  f  r_  f^\T  .  C^V^y  rv„  C^M/.  /'-vA*> 


/>  r  nr  %mp  ^  /•  • '  » P  jk 

j^B^  W  ^  J  ^B/)')^A/y) 


Pa 


,  h'^8 w* 


„™  -  s  J  ^  (*)C  (^>*‘*1  ^yi'-iy 

_^A  q  %  xi  yi 

^mn  2L  tm 

+  £  J /b,  (^)/a,  (^)  ^  j ^B, (y)v^A,  iyfdy 


^mn  ~  £ 


TE  ^ 

%  kZ  ^1^0  '"B,  .’inp 


(t  inp  tmp  J. 

J  ^A^  (^)  dx 


P  '^Up^  Xi 


I  !vAjy)v^. 


Km,=  'Z 


^  flAlt/  Ifw/?  ^ 

j  ^A^  (^)^A^  (^)  ^ 


"PA 
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VIII.5 


Vu 


■A  r  )k  f  * 

^mn  =  ^  ^  W  ^J^A/y)^A/>’)  ^ 

P  9  xi  yi 


TE,  TM-TO-X  MODES  IN  SECTION  A,  TE,  TM-TO-V  MODES  IN  SECTION  B 


Nf 

^mn  ~  ^ 
P 


’^Zb' 


^pb 
nI^x^ 


[C,]  =  [C  +  F  +  Kn 

d+g+l] 

[Cv]  =  [a  +  H-P]'' 

f  ^np  ^mp  ^ 

j  (x)h^  (x)  dx 

yu 


"  f  ,mq  ,mp  ^  r  <<q  <P  * 

+  £  (x)  dx]VB^{y)Vg^{y)  dy 

.  « 


yi 


Hmn=t, 


\tTM  rn  ^  Y 

/Vd  /iK? 

'^2 —  J  ix)VB^  (x)  dx 

P  ‘^^pb 

yu  I'l  iy)i;  i^yf 

^mn  =  X  2  rs,  W  dx]  ' - dy 

P  ^  XI  yi 


D 


ArrA-F|*o£A-*« 


mn 


1 1^- 


^mp 


yu 


P  q  “^O^A 


Nl^ 


\  c  'M  'f^P  ^  f  *P  * 

-  J  (x)Ia^  (x)  dx  j  Vb^  (yWA^  (y)  dy 


Xl 


yi 


Nl^ 


G^„  =  y 

Tnti 

p 


yu 


q  Ka  Xl 


yi 


kh^  ^ 


ii2  * 

r.m 


yu 

f  M/ir  %fnp  a,  f  nT  tp  ^ 

J^B,  (xWa^  ixfdx]lB^(y)lA^iyfdy 


n7 


^mn  ~  X 

P 


P  t"* 

r  Xur^X^A 

N'^Xu 

“  /.  ,nq  ,.mp  *  f  ,q  „p  * 


yu 


q  Xl 
nF  un 


yi 

itmp 


+  X  .,2  J  ^B^  (^)  dx\lB  (y)lAyiy)  dy 


r  k 


h2 

-Upg  Xl 


yi 
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J  dx 


_n7  kl  COE^Ea  X, 

^mn  2^  r 


’’  ’‘‘Kp  'j  .,,  .mp  ,  -q  -p  , 

-  S  ^ - i - Ia,  M'a,  (xfdx  V,(y)Vqiyfdy 

Pa 

q  '^Za  '^za^A  xi  yi 


""  yq  p  p 

^m»  =  X  S  J  *  I  'a,  (>’)%  (>■)  <'>’ 

I  At  jTAAt 

P  Q  '^Za'^Za  xi  yi 


Tr  9 

^^A-Ka  «n 

U  /i  ^A„n  r  "'^P  »k 


VIII.6  TE,  TM-TO-X  MODES  IN  SECTION  A,  TE,  TM-TO-ZMODES  IN  SECTION  B 

[Ci]  =  [C  +  F  +  Kr^[D  +  G  +  L] 

[Cv]  =  [A  +  H-P  +  Sr^  D'*'+G'''  +  L+] 


nI^ 


r  *np  imp  * 

I  (x)  dx 


^  ,2  yu 

D  kK  »2,  r  ''”P  ^  f  "9  'P  * 

+—^  £  W  j  ^y 

^  yi  J 


xf  ty/^o 


/•  anp  amp  ^ 

j  Vb^  (xWb^  (x)  dx 


H=\ 

X^mn  Zu  Ha  TM 

„  im*jn  N„Xu  y„ 

p  Ky  Ky„^c<nq  «mp  *  f  <q  «p  * 

I  W  dAliB,(y)iB,(y)  dy 

P  X,  y,  ] 


^mn  ^ 


wTM  t  2  "  ^  J 


.r/Af  9  W  4  ». 

-  Z  — /.  2 . .  i  W  ^B,(yKB/y)  dy 

r  ^Tb^  ^Tb^  xi  yi 
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A:"  C* 


^  f  "'"P  3k  r  ”P  :k 

S  (xfdxjVg  MVs  M'dy 

4*  1  s,  y, 

A^V  )■. 

+ - h - 2-  J  h,  (x)Ib^  (x)  dx\VB^(y)VB^{y)  dy 

hi  hi  "  ^/  3’; 


klsA-hl  y  ,„p  ^  y^  ,g  ,p  ^ 

^  !—• - W  dxjVgWV  (y)  dy 

N™  q  fflcocj  X,  y, 

^mn  =  X  r  „2  *  1 

P  h^A-^xA^p  „nr  ,mp  ^  ,.r  ,p  ^ 

r  hie  A  .  ^  ^  ^  ^ 


nI^  ,2 


X'  ^QegCOllQ  ^  <nq  ,mp  *  ^,q  ,p  ^ 

L  \  J  ^ J  ^b/^)^a/>')  ‘^y 


^mn  ~  ^ 


xrT’Af  1  2 

^A  q  I^TB„^Za 


nI^  l^n  X 


7,^  r  ^  I*  itr  »p  * 

+  £  -^  J  ^  J ^B,  (>’)^A,  (>’)  ^>' 


-  fcWI  ■' 

'■  ZA 


t^2„  ^u  yu 

TF  ~  2  7S“I'^B  dxUg  iy)I^  (y)  dy 

TE  ^  ^tri  J  °X  J 

9  ^B„  y; 

^mn  ~  ^  TF 

„  W  jt" 

B  '^zd  f  "”'■  "'"/’  *  f  "'■  "P  * 

+  2  —  I  ^B  W^A  W  ^  I^B  (>')^A  W 

(OUq  j  °x  ^x  j  °y  '^y 
.  ’’  yi 


T  2  ’2  *' 

koBA-^xA  ^u 

I  ^mp  J  J.  ,np  ,mp 

TM  - * - J  h,  (x)1a^  (x)  dx 

xtIM  1  tfl  J  X  X 

%  ^ZA  "^O^A  X, 

Qnn  -2^  p  2  *' 

^  N^  k^-hA^p  y  ,^p  ^  yu  ,.^  ^ 

_2  J= - ; - ^J/a  (x)/a  (x)  dx  f  Va  (y)VA  (y)  dy 

hrZn  Kn  J  ^  ^  J  y  y 

9  '^z^  '^Za^A  xi  yi 


nfN^lko^A-kxt  y- 

Fq..=  1  I  ^  lvTAy)V::^xfdxlll(y)lI(yfdy 

,  n  ,  m  J  X  X  J  y  y 


p  q  '^za'^za  Xi 
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VIII.7 


^m«=  I 

P 


ii2 


^0  ^xA 


np 


k^^WflQ 


i*  titip  utnp  ^ 

I  (xWa^  (J^)  dx 


TE,  TM-TO-Z  MODES  IN  SECTION  A,  TE,  TM-TO-XMODES  IN  SECTION  B 
[Ci]  =  [C-F  +  K  +  Mr^[D-G  +  L  +  N] 

[Cv]  =  [A  +  H-Sr^f 


D^-G'^'  +  L'^  +  N'^ 


„p  ,n,p  * 

-  j  Mh,  (x)  dx 


^mn  ,  ff. 


Nf 


^mn  ^ 
P 


2  t*2 


*mp 


kzB  X, 


y  Mrt/7  \mp  ^ 

J  (^)  dx 


+i 


^XB 


mp 


yu 


'^Zb  '^zb  Xi 


X  ,nq  «mp  *  r  '9  "P  * 

\Vb^{x)Vb^  (x)  dxjlB^(y)lBy(y)  dy 
yi 


^mn  ^ 


2  '  ^ 
^0  “  ^xB 


’nq 


~  —  ir^  urn  p 

p  q  ^zr^zb 


Xu  yu 

0  tfia  ttfnp  J.  r  ”P 

J/b^W/b^  U)*&jFB,(>')''B,W  rfv 


Nl" 


Dm„=  1 

P 
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